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Preface 


Although Charles Babbage laid down several ideas for computing “en- 
gines,” the forerunners of today’s computers were largely developed in the 
early 40’s as part of the war effort. The Robinson series cryptoanalytic ma- 
chines developed in England in 1941 spawned many families of computers 
still in use today. The MIT differential analyzer and real-time aircraft simula- 
tion project led to the Whirlwind, and eventually to the immensely successful 
DEC family of PDP computers (Programmed Data Processors). And, of 
course, Eniac, built at the University of Pennsylvania, was the guiding light 
behind Univac, IBM and many other successful manufacturers. 

Many people today poke fun at these early machines and regard them as 
dinosaur-like relics. However, it is interesting to consider that a large-scale 
computer of about 25 to 30 years ago had about the same amount of power as 
a typical personal computer of today. It was generally not as reliable or user- 
friendly as a personal computer, and, of course, cost tens of thousands times 
as much. Why bring this up? 

Because a computer of the 50’s required that the programmer be very 
clever and resourceful to solve problems within the capabilities of the com- 
puter. He did not have vast gobs of memory available, blinding quick calcula- 
tion speed, or random disk access. In other words, he had about the same 
problem to face as you do with your personal computer. 

I do not mean to imply that your personal computer is not a full-fledged 
computer. It certainly is just as much a computer as a room-filling giant of to- 
day. However, because of the relatively small memory, it cannot store a large 
data base. Nor is it suitable for extensive word processing or massive calcula- 
tions. Highly detailed graphics are best left to other machines as well. 

What can we learn from the computing pioneers of the 50’s that will help 
us today? Perhaps most important is the discipline of thoroughly analyzing a 
problem, breaking it down into manageable steps, and solving it a step at a 
time. It is also important to determine what can be done “off line” and what 
must be done on the computer. 
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That is what this book is all about. While it has more than 50 ready-to-run 
programs, the main thing you should look for from the book is an approach 
to solving problems — big and small. Some of the problems demonstrate the 
capability of the computer; others identify its shortcomings. It is important to 
be familiar with both the strengths and weaknesses of your tools so you can 
recognize the types of jobs for which they are suitable (and not suitable). 

This book focuses primarily on mathematical and educational applications 
for the computer. There are many other excellent sources of information 
about other applications and for making best use of your personal computer. 
Using the approaches described in this book should enable you to easily con- 
vert programs and use applications from other books and magazines such as 
Creative Computing. 

This book is designed to be read with a working computer at hand. While 
there is textual material to be read, the most important things are the experi- 
ments and problems to be tried with your own computer. The book raises 
many questions for which you should try to find answers. There are no an- 
swers to these questions and problems in the back of the book; you should be 
able to discover the answers as you work the problems out on your computer. 

You will be able to incorporate many of the routines and approaches in the 
book into programs of your own as you use your computer to deal with “real 
world” problems. Other programs will simply point you in the right direc- 
tion. And some of the programs in the book are just plain fun. Learn. Experi- 
ment. Have fun! 

Morristown, New Jersey 

March 1983 David H. Ahl 
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and Practice 


Throughout life, there are certain things that simply must be memorized. 
Obvious things that fall into this category are the addition and multiplication 
tables, the spelling and meaning of words, how to tell time, and the monetary 
system. 

However, depending upon one’s chosen profession, there are many other 
things to be memorized. A doctor must know what diseases match what 
symptoms. A chemist must know the gas laws, the properties of elements and 
so on. A pilot must instantaneously know the meaning of readings on scores 
of instruments. 

To memorize a set of facts, you must go over them again and again and 
keep trying different variations. Here is where the computer comes in. It is 
able to present randomly scores of different problems to you for as long as 
you wish. Some programs will automatically adjust to your level of com- 
petence and will grade you; other programs simply present the problems and 
leave the grading up to you. 

There are four programs in this chapter and two in the Science chapter 
which present material in a drill and practice format. Examine the methods 
used in these programs and then make up some drill programs of your own 
for subjects with which you are having trouble, or make up programs for 
other members of your family. 


Addition Practice 


This program demonstrates a simplified addition drill and practice routine. 
This type of drill is sometimes called computer assisted instruction (CAI), al- 
though CAI can also apply to tutorial and other approaches as well. 

When the program is run, it will first ask “No. of digits?” You enter a 
number and each addend will contain that many or fewer digits. 

The program will present any number of practice problems that you spec- 
ify. The program presents each problem in turn. The program will not pro- 
ceed to the next problem until the current one has been answered correctly. 
After the last problem is answered correctly, the score is printed with an 
appropriate comment. 

There are many improvements and extensions possible in a program like 
this one. For example, you might want to modify the program so it tells the 
user the correct answer after a problem has been answered incorrectly two (or 
three) times. 

A more complicated modification would be to change the program to 
present different kinds of arithmetic problems such as subtraction, multiplica- 
tion, and division. 

Some of these modifications have been made in the next program. 


5 POR I«1 TO VAL (RIGHTS (TIMES .2)) 

10 DUMMY*RND( 1 ) : NEXT I 

20 CLS: PRINT TAB(11 ) i "Addition Practice” : PRINT 

30 INPUT "Number of digits" ;N 

40 INPUT "Number of problems" ;P 

50 W-0 

60 X=X+1 

70 Q=*0 

80 A=INT( 10“N*RND( 1 ) )+l 
90 B=INT(10“N*RND(1))+1 
100 CLS 

110 PRINT TAB( 8-INT( LOG(A)/LOG( 10 )+l ) ) ; A 
120 PRINT TAB(7- INT( LOG( B )/LOG( 10 )+l ) ) j "+" ; B 

130 PRINT " " 

140 c*A+B 

150 PRINT TAB(7-INT(LOG(C)/LOG(10)+1)) ; 

160 INPUT G : PRINT 
170 IP G-C THEN 240 
180 Q-Q+l 

190 IP Q>1 THEN 210 
200 W-W+l 

210 PRINT "What? Try again... "i 

220 GOSUB 340 

230 COTO 100 

240 PRINT "Right! "s 

250 IP X>-P THEN 290 

260 PRINT "Here's another..." 

270 GOSUB 340 
280 GOTO 60 

290 CLS: PRINT "You got" ;P-W;"correct the first time." 
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3*0 IP W< . 22*P THEN 33* 

31* PRINT "But you missed" sW 
320 COTO 35* 

33* PRINT "Good work!" 

34* FOR Z-l TO 999 .NEXT Z: RETURN 
350 END 


29 
+ 488 


517 

Right! Here’s another... 


654 
+ 451 


1005 

What? Try again... 


654 
+ 451 


1105 

Right! Here's another... 


826 
+ 816 


1642 

Right! Here’s another... 


582 
+ 37 


619 

Right! Here’s another... 


254 
+ 369 

613 

What? Try again. . . 


254 
+ 369 


623 

Right! 

You got 3 correct the first time. 
But you missed 2 
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Addition Practice, 
Adjusted by Grade Level 


One of the major disadvantages with many drill and practice exercises is 
that they tend to be either boring or frustrating, depending upon the ability of 
the user relative to the level of the material. To compensate for this, a method 
is needed which will adjust the difficulty of the problems to the ability of the 
user. 

Ideally, such a system would weigh the most recent performance most 
heavily but would not ignore previous performance. It should allow a user to 
advance to more difficult problems than his current mastery level. It should 
also continue to give some practice on problems already mastered. 

Some commercial software packages approach these goals along traditional 
lines, i.e., determine in which type of problems a student should receive prac- 
tice by using a complicated computer managed instruction score recording 
and adjustment system. 

The approach here is more innovative; it uses a single measure for each 
type of problem — call it “estimated grade level” — which meets all of the 
objectives stated above. 

How does it work? The most recent problem presented counts 10% of the 
overall score if it was answered correctly and was over the current user grade 
level, or if it was answered incorrectly and was under the current user grade 
level. Otherwise it is ignored. This may be easier to visualize in the form of a 
chart: 

Answer 

Right Wrong 


Problem 


Higher than 
grade level 

Lower than 
grade level 


Raise student 
grade level 

Ignore 


Ignore 

Lower student 
grade level 


At first glance this might look complex and somewhat goofy, however, 
what it really means is that a student is rewarded for doing a problem beyond 
his grade level but he is not penalized if he cannot do it. On the other hand he 
is penalized if he cannot do a problem lower than his grade level, but is not 
rewarded for doing one lower. 

Each problem affects the estimated grade level a little bit, with the most 
recent problems being weighed the most heavily. If the current grade level of 
a student is L and the level of the most recent problem to be averaged in is P, 
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then the averaging formula is simply: 

L = 9L + .IP 

The remaining task before a program can be written is to assign a grade 
level to each problem presented. Unfortunately, this will vary depending 
upon the local school system, the textbook used, and the teaching method. 
Also a huge data base can not be stored in a small computer, so it is desirable 
to devise a simple method of determining grade level for different problems. 
One straightforward approach is to present problems up to one-half a grade 
level over and under where the student currently is. Thus the overall range of 
problems for a student at grade level 3.2 would be 2.7 to 3.7. 

How do we generate the right problems? Consider one type of skill, vertical 
addition. It is normally introduced in the first grade and continues through 
Grade 4 (actually 4.9). The simplest problem in this program is 1 + 1 and the 
most difficult is 999 4- 999. Since learning is not a linear process (it is slow at 
first, and then progresses rapidly), an exponential formula can be used. For 
example: 

Addend = 1.73 x (Grade level) 4 
or 


Grade level = V Addend/1.73 

This means that students at various grade levels will be working with the 
following maximum addends: 


Grade level 

Addend 

1.0 

1 

2.0 

27 

3.0 

140 

4.0 

442 

4.9 

997 


Now it is a relatively straightforward, although somewhat tedious, matter 
to tie all these elements together in a computer program. 

A few notes about the program. The variable G2 is the problem grade level 
that is always within one-half of a grade level of the current student level, Gl. 
The complicated mess in Statement 310 produces a moving grade level 
average. 

The recording of the grade level and carrying it over to the next lesson is a 
manual process. On a computer system with a permanent mass storage de- 
vice, this would be kept on the system. 

There are many possible changes and extensions to this program. For 
example it could present different types of problems such as horizontal addi- 
tion, vertical and horizontal subtraction, as well as multiplication, division 
and fraction problems. 
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5 FOR 1=1 TO VAL (RIGHTS (TIMES ,2) ) 

10 DUMMY =*RND ( 1 ) ; NEXT I 
20 W=0 

40 tnph? R »!! T ” Hi ' T ° St0P ' enter 9999 as answer." 

40 INPUT Your current grade level";Gl 

50 G2=G1- . 5*RND( 1 ) 

60 R=INT(2+1.73*G2~4) 

70 A*INT( 100*G2*RND( 1 ) ) 

80 IF A>R THEN 70 

90 B=R-A 

100 PRINT TAB ( 8-INT( LOG( B )/LOG( 10 )+l ) ) ; B 

120 PRINT ^ AB ^:^ NT(L ° G(A)/LOG ( 10 ) + l>)i"+ ,, iA 
130 R*A+B 

15» iSp 5? ^-INTO.OCCRiAOm.Hl,,, 123 

160 IP G=9999 THEN 340 + 119 

170 IF G=R THEN 280 IT!" 

180 W=W+1 242 

190 IP W>1 THEN 220 Correct! 

200 print "Wrong, try again..." *. w 

210 GOTO 100 Here's another... 

220 PRINT "You missed that one twice." 

230 PRINT "Correct answer is":R . Q 

240 GOSUB 380 *? 

250 W=0 + 111 

260 IF G2>G1 THEN 320 ”""~ 

270 GOTO 310 r . rr Jr 

280 W*0 Correct! 

290 PRINT "Correct! ": GOSUB 380 Here's ennfh.r 

300 IF G2<G1 THEN 320 Here s another... 

310 Gl= . 9*G1+ . l*SQR(SQR(R/2/l . 73) ) 

320 PRINT .PRINT "Here's another..." 

330 GOTO 50 

340 CLS + 163 

350 PRINT "Okay. So long for now." Ill' 

360 G1“(INT(100*G1))/100 ,, 

370 PRINT "new grade level is";Gl:END ror >g, try again.. 

380 FOR Z-l TO 900: NEXT Z. -RETURN + ^ 


49 

+ 111 
160 

Correct! 

Here’s another. . . 


67 
+ 163 


Hi. To stop, enter 9999 as answer. 
Your current grade level? 3.5 


You missed that one twice 
Correct answer is 230 

Here's another. . . 


193 
+ 50 


243 

Correct! 

Here’s another. 


61 
+ 171 

232 

Correct ! 

Here’s another 


230 
+ 27 


106 
+ 58 


257 

Correct! 


164 

Correct ! 


Here ' s another . . . 


Here’s another. 
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75 
+ 67 


50 
+ 146 


9999 

A Okay. So long for now. 
New grade level is 3.43 


Hi. To stop, enter 9999 as answer. 
Your current grade level? 3.5 

47 

+ 199 


246 

Correct! 

Here ' s another . . . 

Ill 
+ 57 


168 

Correct ! 

Here's another . . . 

55 
+ 155 


210 

Correct! 

Here's another . . . 
5 

+ 153 


158 

Correct! 

Here's another. . . 

191 
+ 44 


235 

Cortect! 

Here ' s another . . . 

187 
+ 9 


196 

Correct! 

Here ' s another . . . 


196 

Correct! 

Here’s another. . . 

148 
+ 11 


159 

Correct ! 

Here ' s another . . . 

101 
+ 124 


225 

Correct! 

Here ' s another . . . 

161 
+ 75 


226 

Wrong, try again... 
161 
+ 75 


246 

You missed that one twice. 
Correct answer is 236 

Here’s another . . . 

Ill 
+ 40 


151 

Correct! 

Here’s another. . . 

64 
+ 144 


9999 


Okay. So long for now. 
New grade level is 3.43 
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Time/Speed/Distance Problems 

As well as being able to present simple numerical drill and practice prob- 
lems, the computer can present word problems for solution as well. 

In this program the formula relating time, speed and distance was applied 
to a problem involving both a car and train. 

The problem can be stated as follows. A car traveling C miles per hour 
(computer generates an integer 40 through 65) can make a certain trip in D 
hours (computer generates an integer 5 through 20) less than a train traveling 
at T mph (computer generates an integer 20 through 39). How long does the 
trip take by car? When the two simultaneous equations are solved they pro- 
duce the single equation for the answer shown in Line 90. 

Notice the calculation in Line 100. This calculates the percent difference 
between the actual answer and the one entered by the user. 

Notice also that the computer calculates the correct answer in Line 90 and 
prints it (on the screen) in Line 180. This answer may have many decimal 
places; as the program is written it is rounded off to two decimal places. If the 
0.5 was not added in Line 100, the number would be truncated and not 
rounded off. This is an important calculation and one which you will find in 
many other programs throughout the book. 

Consider other problems that can be used as the basis for this kind of drill 
and practice exercise. Teachers, for example, might wish to have students 
write drill and practice programs on their own. Different problems could be 
given to one or a small group of students to serve as the basis for a program. 

After the programs are written, students can try out the programs of other 
class members. This approach ensures that each student not only understands 
the type of problem assigned to him, but also gets practice in solving other 
problem types as well. This is an effective technique for stimulating interest as 
well as for learning how to solve word problems. 
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5 FOR I«1 TO VAL(RICHT$<TIME$,2)) 

10 DUMMY“RND( 1 ) j NEXT I 
20 C-INT(RND(l)*25)+40 
30 D«INT(RND(l)*15)4-5 
40 T«INT(RND(l)*19)+20 
45 CLS 

50 PRINT "A car traveling at”jCs*'mph makes a" 

60 PRINT "trip in" iDj "hours less than a train" 

70 PRINT "traveling at" ; T,"mph. How long is" 

80 INPUT "the car trip in hours” : A 
90 V«D*T/(C-T) 

100 £-INT(ABS((V-A)*100/A)+.5) 

110 IF E>5 THEN 160 

130 PRINT "Good. Answer is within" ;E S "percent." 

140 GOTO 180 

160 PRINT "Sorry. You were off by" ;E; "percent." 

180 PRINT "Correct time is" s INT(100*V)/100 , "hours." 
190 PRINT: INPUT "Try another (Y or N)" f z$ 

200 IF Z$*"y" OR Z$-”Y" THEN 20 
210 PRINT "Okay. So long for now." 

220 END 


A car traveling at 42 mph makes a 
trip in 16 hours less than a train 
traveling at 30 mph. How long is 
the car trip in hours? 48 
Sorry. You were off by 17 percent. 
Correct time is 40 hours. 

Try another (Y or N)? 


A car traveling at 63 mph makes a 
trip in 14 hours less than a train 
traveling at 28 mph. How long is 
the car trip in hours? 11 
Good. Answer is within 2 percent. 
Correct time is 11.2 hours. 

Try another (Y or N)? n 
Okay. So long for now. 

Ok 
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Kinematics Problems 


Kinematics relates to the dynamics of motion of bodies apart from consid- 
erations of mass and force. Like many other types of problems, these can be 
generated and presented by the computer to provide practice in solving them. 

This program presents a simple kinematics problem for solution. The com- 
puter generates a new value for each problem. The problem is as follows. A 
ball (or any other object) is thrown up at velocity V meters per second (com- 
puter generates an integer between 5 and 40). The user must then calculate 
three factors about the resulting flight of the ball: maximum height, time until 
it returns, and velocity after T seconds (computer generates a time less than 
the total flight time). 

The key benefit in using a computer to present problems of this type is 
motivation. The calculations required of the user are no different than those 
in the back of a chapter in a book or those on homework assignments. How- 
ever, answering them when they are presented by the computer seems to 
make it more of a challenge and, frankly, more fun. 

The computer program checks each of your responses to see if it is within 
15% of the correct answer. If it is, your answer is considered correct. You 
may wish to change this percentage to require more or less accurate calcula- 
tions. You may also wish to change the computer calculations to round off to 
one or two decimal places. 






5 FOR I“1 TO VAL(RIGHT$(TIME$,2) ) 
10 DUMMY*RND ( 1 ) : NEXT I 


20 Q“0 

30 V=6+INT(35*RND(1)) tv „ 

40 PRINT ”A ball is thrown upwards at a velocity 

50 PRINT "of" 5 Vi”meters per second.” 

60 A* 05*V~2 

70 INPUT "How high will it go”;G 
80 GOSUB 250 

100 A INPUT "How many seconds until it comes down ;G 
110 GOSUB 250 

130 T=2+INT(2*V*RND(l))/10l 
140 IF T>V/5 THEN 130 

150 A“V-10*T if 

160 PRINT "What is its velocity at";Tj 

170 INPUT "seconds" ;G 


190 GOSUB 250 _ , - „ 

200 PRINT: PRINT Q fright out of 3. 

210 IF Q>1 then PRINT "Not bad. 

220 PRINT: INPUT "Try another (Y or N) jss 
010 IP ZB*"v" OR Z$*"Y" THEN 20 


240 END __ _ 

250 IP ABS C(G-A)/A)<-15 THEN 280 
260 PRINT "Not even close..." 

270 GOTO 300 

280 PRINT "Pretty good." 


290 Q-Q+l „ 

300 PRINT "Correct answer is ;A 


310 RETURN 


A Ball Is thrown upwards at a velocity 
of 40 aeters per second. 

How high will It go? 80 
Pretty good. 

Correct answer Is 80 

How nany seconds until It coses down? 8 
Pretty good. 

Correct answer Is 8 

What is Its velocity at 2.8 seconds? 10 
Not even cl ose. . . 

Correct answer Is 12 

2 right out of 3. 

Not bad. 

Try another (Y or N)? N 
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Problem Solving 


In many courses in school, the textbooks present a great variety of devices 
for solving the problems that have been neatly grouped together at the end of 
each chapter. Typically these devices consist of formulae, equations, rules and 
theorems. After a careful study of these devices, teachers give exams which 
test your ability to recall them. 

But what do you do if you are faced with the more realistic situation of not 
being told what device is likely to solve which problem or, worse yet, of hav- 
ing forgotten how to use a technique altogether? Is all hope lost? Of course 
not, although some people seem to believe that it is. 

In this chapter, several of these nasty devices mentioned above are pre- 
sented. For example, there are devices written into computer programs that 
will solve a quadratic or exponential equation and others that will calculate 
the roots and draw a plot of any function. 

However, you must remember that while these devices are useful in solving 
certain problems, the really important thing is to understand the underlying 
logic and approach. Then when it comes time to solve real world problems 
you will be better prepared to face them. Incidentally, many of the methods 
and devices presented in this chapter are used in later chapters to solve other 
kinds of problems. 
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How Many Tickets? 


Here is a problem. At a school raffle to raise money, the organizers have as 
a prize an electronic game for which they paid $18.00. To add interest to the 
raffle, the organizers have decided to sell tickets for an amount (in cents) 
equal to the number on the ticket for tickets numbered 1 to 50. For ticket 
numbers over 50, the price is 25 cents each. The organizers want to know 
how many tickets they must sell to exactly break even. 

It is probably easiest to visualize a problem of this sort with a flowchart. In 
the flowchart, T will equal the total money collected and will increase as 
more tickets are sold. The ticket number is N. When T equals or exceeds 
$18.00, N will be the answer. 

Note that the flowchart has two logical branching points (if statements in 
the program). The first compares the current ticket number to 50; if it is less, 
the ticket number is added to the total whereas if it is greater than 50, the 
total is increased by 25 cents. 
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The second branch point compares the total amount collected, T, to $18.00 
(actually 1 800 cents). If T is equal to or greater than 1 800, the break even 
point has been reached and the values of N and T are printed (on the screen). 

A problem of this kind can be done by hand, however, because of the repet- 
itive additions it is quite tedious. Also, doing it by hand frequently leads to an 
answer of 72 rather than the correct answer of 71. Try it yourself and see 
what you get. 


10 T = 0 

20 FOR N=1 TO 2000 
30 IF N< = 50 THEN 60 
40 T=T+25 
50 GOTO 70 
60 T=T+N 

70 IF T>=1800 THEN 90 
60 NEXT N 

90 PRINT N; "Tickets sold" 

100 PRINT "Collected $* ; T/100 


71 Tickets sold 
Collected $ 18 


Many problems can be solved quickly and correctly with a computer using 
logical analysis and a flowchart. More complex problems may have to be bro- 
ken down into additional steps and require a longer flowchart, but the ap- 
proach is fundamentally the same. 

Here are two problems for you to solve. 

The diameter of a long-playing record is 12 inches. The unused center has a 
diameter of 4 inches and there is a smooth outer edge 1/2 inch wide around 
the recording. If there are 91 grooves to the inch, how far does the needle 
move during the actual playing of the recording? 

A movie theater charges $2.50 for an adult admission and $1.00 for a child. 
At closing, the cashier counted 385 ticket stubs and had $626.50 in cash. 
How many children entered? 
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Drinking and High Blood Pressure 

This program illustrates how several simple equations can be put in a com- 
puter program to solve a more difficult overall problem. 

Here is the problem. In a survey of 1000 adults, it was found that 35 had 
high blood pressure. Of those with high blood pressure, 80% drink 15 oz. or 
more of alcohol per week. Of those without high blood pressure, 60% drink a 
similar amount. What percent of drinkers and non-drinkers have high blood 
pressure? 

This problem requires the solution of several simple equations. They could 
all be combined into one large equation, but it may be easier to understand 
the approach (and change variables later on) by writing a program with five 
separate equations. 

If H equals the number of people with high blood pressure, then H = 35 
(Line 10). Letting HI equal the number of people with high blood pressure 
who drink leads to HI = .8 x H (Line 20). 

Letting LI equal the number of people with low blood pressure who drink 
yields LI = .6 x (1000 - H). Then, the total number of drinkers, D = HI + 
LI. 

Finally, the percentage of drinkers with high blood pressure is X = HI x 
100 /D. 

The program solves the problem in a jiffy. The solution for this type of 
problem can be easily written directly in Basic without any need for a flow- 
chart or detailed analysis. Recognizing this type of problem readily will save 
a great deal of pencil pushing time. 


10 H=35 

20 HI* . 8*H 

30 Ll=. 6*(1000-H) 

40 D=H+L1 
50 X=H*100/D 
60 X1*H*100/1000 

70 PRINT ■ High Blood Pressure" 
80 PRINT "Population *;X1}"%" 

90 PRINT "Drinkers * ; X } "%* 


High Blood Pressure 
Population 3.5 % 
Drinkers 5.70033 % 
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Here is a problem that doesn’t require a single equation but makes use of 
the approaches discussed so far in this chapter. Can you solve it with three 
Basic statements? 

In early January, a shopkeeper marked down some calendars from $2.00 to 
a lower price. He sold his entire stock in one day for $603.77. How many did 
he have? 

Here’s another easy one. A town in India has a population of 20,000 peo- 
ple. Five percent of them are one-legged and half of the others go barefoot. 
How many sandals are worn in the town? 
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Two Simultaneous Equations 


So far in this chapter, only problems with linear equations have been 
considered. But the computer can be used to solve much more difficult equa- 
tions. In fact, it is in problems involving second and third degree equatio , 
exponentials, and the like where the computer really starts to pay off. Con- 
sider the following two simultaneous equations: 


2X = 16y 3 X = 27y 

It is not at all easy to solve these two equations by hand. But a simple Basic 
program can be written to solve the equations using trial and error. This is 
sometimes referred to as a brute force approach because every possible 
combination of numbers between an upper and lower limit is tried until a 
solution is reached or until the program runs out of values. 


10 PRINT "Solves two simultaneous equations' 
20 FOR X=1 TO 10 

40 i? R 2‘'X<>16*Y/3 OR 3 a X<> 27*Y THEN 70 
50 PRINT * x* " ; X, "y="{Y 
60 END 
70 NEXT Y 

80 NEXT X , _ „ 

90 PRINT "No Integer solution" 


Solves two simultaneous equations 
X* 4 y= 3 


In this particular case, a solution is reached rather quickly with x = 4 and 
y = 3. However, if in the second equation the y coefficient is changed slightly 
from 27 to 28, the computer will try 10,000 possible solutions before finally 
concluding that no integer solution exists— at least within the range of 0 to 
100. Warning: this will run for a very long time. 


40 IF 2 A X<>16*Y/3 OR 3 A X<>28*Y THEN 70 


Solves two simultaneous equations 
No Integer solution 
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Although the trial and error (brute force) approach is widely used, it is 
highly inefficient. In general, a systematic or guided trial and error approac 
is preferable to one that simply tries every possible solution. However, for 
some problems the simple “try every value” approach may be appropriate. (A 
comprehensive discussion of trial and error approaches can be found on pp 
36-40 of Computers in Mathematics: A Sourcebook of Ideas.) 

Three problem solving approaches have been discussed so far. Remember- 
ing them, how would you do this problem? A boy and his sister visited a farm 
where they saw a pen filled with pigs and chickens. When they returned 
home, the boy observed that there were 18 animals in all, and his sister re- 
ported that she had counted a total of 50 legs. How many pigs were there in 

the pen? 
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Quadratic Equation Solver 


For any values A, B, and C of a first degree quadratic equation (Ax 2 + Bx 
+ C = 0), this program will compute the roots of the equation. The solution 

is based on the quadratic theorem which solves for roots with the following 
formula: 


-B±V / B 2 -4AC 

2A 


Assuming that A, B, and C are real numbers, the following principles 
apply: 

1. If B 2 - 4AC is positive, then the roots are real and unequal. 

2. If B 2 - 4AC equals 0, then the roots are real and equal. 

3. If B 2 - 4AC is negative, then the roots are imaginary and unequal. 

The program takes into account all these possibilities and correctly identi- 
fies the type of roots along with their values for any set of coefficients. 

Is this program useful by itself? Except for solving quadratic equations for 
algebra class, probably not. However, as a routine in a larger program to 
solve quadratic equations that might be encountered, it could be very useful. 


10 PRINT 
20 INPUT 
30 R=8*B- 
40 IF A<> 
50 PRINT 
60 PRINT 
70 IF R<0 
60 IF R=0 
90 PRINT 
100 GOTO 
110 PRINT 
120 PRINT 
130 PRINT 
140 PRINT 
150 INPUT 
160 IF A$ 
170 END 


"Quadratic Equation Solver" 
A,B,C";A,B,C 
4*A*C 
0 THEN 60 

"First degree equatl on. " : GOTO 150 
"Roots are 
THEN 120 
THEN 110 

(-B+SQR(R) )/2*A ; (-B-SQR(R) )/2*A 
150 

-B/2*A : GOTO 150 
" 1 iagl nary: * 

(-B/2*A) ; " + ";SQR(-R)/2*A; "*t « 
(~B/2*A) j "-";SQR(-R)/2*A; **1 ■ 
"Another (Y or N)" ; A$ 
s "y" OR A$="Y" THEN 10 


Quadratic Equation Solver 
A,B,C ? 1 , 2 , 1 
Roots are -1 
Another (Y or N)?y 


Quadratic Equation Solver 
A, B, C ? 0 , 2 , 4 
First degree equation. 
Another (Y or N)?y 


Quadratic Equation Solver 
A, B, C ? 2 , 5 , 2 
Roots are -2 -8 
Another (Y or N)?y 

Quadratic Equation Solver 
A,B,C ? 4 , 2 , 4 
Roots are liaglnary: 

-4 + 15.4919 *1 
-4 - 15.4919 *1 
Another (Y or N)?n 
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Exponential Equation Solver 


Another general routine for solving a particular type of equation is this one 
to solve for an exponent in an exponential equation. 

Given the values A, B, m, and n, this program will solve for x in any 
exponential equation of the form: 

A mx+n = B 

For example, the program will solve any of the following problems: 

1. 5 X = 40 

2. 5 3X+I = 7.6 

3. 17 x * 3 = 8.12 

4. 11>- 2X = 247 

If you were to solve an exponential equation by hand, you would probably 
go through the following steps: 

5 X = 40 

log 5 X = log 40 

x log 5 = log 40 

x = log 40/log 5 

x = 1.6021/.6990 = 2.292 

However, in more generalized form, the solution for x is: 

log B 

X = - lo K A 

M-(N/M) 

Note that the program to solve this problem is actually divided into two 
sub-programs. The first is a data loader program. It requires that data be en- 
tered in the following order: A, B, m, and n for each equation to be solved. If 
a coefficient is not present, it must be entered as a zero. 

The data for the four problems listed above were entered into the four data 
Statements (120-150). 

The program as it is presented here can be improved in several ways. First, 
it always solves the same four equations. How can you generalize it to solve 
for other equations? Second, if you were to make use of this routine in an- 
other program, you would probably not be able to use a read statement; how 
could you get rid of it?. 
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10 CLS: PRINT"Exponent1al equation solver" 

20 PRINT : PRINT • A B * n 

30 READ A, B, M, N X 

40 X=(LOG(B)/LOG(A))/M-(N/M) 

50 X=INT(. 5+1000*X ) / 1 000 

70 Q^5 T 30 ;TAB(6) ' BjTAB(14)iMiTAB(22)|N!TAB(30)iX 

120 DATA 5,40,1,0 

130 DATA 5, 7. 6, 3,1 

140 DATA 17,8.12,1,-3 

150 DATA 11,247,-2,1 


Exponential equation solver 



You won’t believe this, but Zarg can actually add without a computer. 





Roots of Any Function 

This program will find the roots of a function, any function! The function 
may be linear, quadratic, cubic, trigonometric or any combination as long at 
it can be represented in the Basic language. The program as it appears here 
finds the roots between -20 and 20 although you can change these boundaries 

in Statement 40. 

The method used involves evaluating the function at small incremental 
intervals, finding places where the value of the function changes sign and 
then, by successive approximations, finding the zero point. This approach 
borrows from Newton’s method in the final narrowing down but, unlike 
Newton’s method, will not fail to converge in the event one makes an unlucky 
first guess. 

Before running the program you must first type in your function in state- 
ment 90. For example, 

DEF FNA (X) = 2*XT3+11*Xt2-31*X-180 

DEF FNA (X) = X- 4 

DEF FNA (X) = SIN(X) - . 5 

You may have to refer to the Basic manual with your system to see exactly 
how a function should be stated. 

The routine used in this program is very powerful and could possibly be 
used as a subroutine in many other programs. 

How f- fln you use this program to help you solve this problem for x? 


x = 
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10 PRINT "Roots of any function. Define" 
20 PRINT "function In subroutine 260." 

30 PRINT: Zl=-200 

40 FOR I=-19. 9 TO 20 

50 X = I : GOSUB 260: T1=FNX 

60 X=I+1 : GOSUB 260: T2=FNX 

70 IF SGN(T1)=SGN(T2) THEN 240 

80 K= I : 0= I + 1 

90 IF T1<T2 THEN 130 

100 Z=I 

110 K = I + 1 

120 J=Z 

130 Z*(K+0) 12 

140 X=Z: GOSUB 260: T3=FNX 

150 IF T3<0 THEN 180 

160 0=Z 

170 GOTO 190 

180 K = Z 

190 IF ABS(T3 ) > 5E-05 THEN 130 
200 Z=SGN(Z)*INT(ABS(Z)*10000+.5)/10000 
210 IF Z=Z1 THEN 240 
220 PRINT "f (";Z* " ) = 0 " 

230 Z1=Z 
240 NEXT I 

250 PRINT "Prograi f 1 nl shed" : END 
260 FNX=SIN(X)-. 5 
270 RETURN 


260 FNX=SIN(X)-.5 


260 FNX=X-4 


Roots of any function. Define 
function In subroutine 260. 

f (-18. 3259 ) *0 
f (-16. 2315 )=0 
f (-12. 0428 )=0 
f (-9. 9483 ) *0 
f (-5. 7596 ) *0 
f (-3. 6651 ) =0 
f( .5236 )=0 
f( 2.6179 ) *0 
f( 6.8068 )=0 
f( 8.9011 ) =0 
f( 13.09 ) *8 
f( 15.1644 )=0 
f( 19.3732 )=0 
Progra* finished 


Roots of any function. Define 
function In subroutine 260. 

f( 4 ) = 0 

Prograi finished 


260 FNX=2*X"3+11*X"2-31*X-180 


Roots of any function. Define 
function In subroutine 260. 

f (-5 ) =0 
f (-4. 5 ) s 0 
f( 4 )=0 

Prograi finished 
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Plot Any Function 


Here is a nifty program that will produce a plot of any function on the 
screen or printer. 

Before running this program, you must type in your function in line 120. 
Like the previous program, which finds the roots of any function, this one 
will plot any function. You must tell the program between what values you 
want the function plotted, i.e., a minimum and maximum value of the x co- 
ordinate. You also input the x plotting increment you wish. 

As it appears here, the program does not allow the user to select the y co- 
ordinates; the program plots y values between -30 and 30. 

Here are several functions you might want to try plotting: 


Function 

X Limits 

Y Limits 

Y = SIN(X) 

-6 . 3 

6 . 3 

-1.1 

1 . 

Y = SIN(X) 

0 

314 

- 1 

1 

Y = EXP ( S 1 N ( X ) ) 

-6 . 3 

6 . 3 

0 

3 

Y = 910/X) * S 1 N ( X ) 

2 

24 

-2 . 5 

5 

Y = 3 * S 1 N ( X ) +X 

-12.6 

12.6 

- 15 

15 

Y = Xf 2- 1 5 

-6 

6 

-20 

30 

Y = XT2+X+20 

-7 

8 

-38 

22 

Y = 30*EXP( -X*X/ 100) 

-25 

25 

0 

30 


The last function listed is the one plotted in the sample run with the pro- 
gram. Exponential functions are a great deal of fun and sometimes lead to 
unexpected and interesting results, particularly when combined with trigo- 
nometric functions. Experiment! Have fun! 


10 CLS:PRINT "Plot any function. Define" 
20 PRINT "function In Line 120." 

30 INPUT "Lower & upper X limits"? XL , XU 

40 INPUT "Lower & upper Y 1 1 ml ts" ; YL, YU 

50 DX=XU-XL 

60 SP=DX/240 

70 DY=YU-YL 

80 SY=DY/64 

90 CLS: SCREEN 0,0 

100 FOR X=XL TO XU STEP SP 

110 XP= ( X-XL ) /SP 

120 Y=3*SIN(X)+X 

130 YP = 64- ( ( Y-YL ) / S Y ) 

140 IF YP<0 THEN YP=0 
150 PSET (XP.YP) 

160 NEXT X 

170 IF SQN( XU) =SGN( XL ) THEN 220 
180 XA=-XL/SP 
190 FOR Y=0 TO 63 
200 PSET ( XA, Y ) 

210 NEXT Y 
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220 IF SGN( YU ) =SGN( YL ) THEN 270 
230 YA=64+YL/SY 
240 FOR X=0 TO 239 
250 PSET ( X, YA) 

260 NEXT X 
270 GOTO 270 
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Sets and Repetitive Trials 


For solving relatively simple problems, the computer may not be any help 
at all. In fact, it may take more time to write a program to solve a problem 
than it would to solve it by hand or with a calculator. This chapter should 
help you recognize problems that are suitable for computer solution and 
those that are not. 

In some problems, you may think that the computer will not be any help. 
However, one thing writing a computer program will always do is force you 
to reason out the approach to solving the problem logically and precisely. The 
computer can’t solve problems unless it is told exactly how to proceed; hence 
you must understand a problem completely before you can program it for the 
computer. 

Several of the sections in this chapter discuss sets of data. While the sets 
used in the examples have relatively few elements or values, you should bear 
in mind that real world problems often have thousands or millions of pieces 
of data, and the only practical way to solve problems of this size is with a 
computer. For example, consider how you would most efficiently schedule 
the shipments on a railroad train leaving Boston with 4000 diverse cargoes 
bound for Phoenix and 780 points in between. Now consider that there are 
200 freight trains per day leaving Boston. Now add to that the 10,000 other 
trains leaving other cities every day that must use the same network of track 
and you can see that dealing with real world sets of data is no easy task. 
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Group of Girls and Boys 

In the previous chapter we said that there may be some problems for which 
brute force trial and error is appropriate. This would be the case if the prob- 
lems were relatively small and trying every possible solution would not tie up 
a great deal of valuable computer time. Here is a problem involving two lin- 
ear equations that lends itself to a trial and error approach. 

The problem is as follows: When 1 5 girls leave a group of boys and girls, 
there are two boys for every girl (lucky girls). Next, 45 boys decide to leave; 
then there are 5 girls for every boy (lucky boys!). How many girls were there 
in the group before anyone left? 

Before rushing to the computer, you must recognize that this problem re- 
quires the solution of two simultaneous equations. If G equals the original 
number of girls and B the original number of boys, then the two equations 
are: 


(G - 15) x 2 = B 
(B - 45) x 5 = (G - 15) 


10 PRINT "Girls and boys problem' 1 

20 PRINT "Program is S-L-O-W. Be patient." 

30 1 = 0 

31 Gl=15 

32 81=45 

40 FOR G=1 TO 100 
50 FOR 8=1 TO 100 
60 1=1+1 

70 IF 2* ( G-Gl ) <>B THEN 110 
80 IF 5* ( B-Bl ) < > ( G-Gl ) THEN 110 
90 PRINT Gj “Girl s M ; B; "Boys" 

100 GOTO 140 
110 NEXT B 
120 NEXT G 

130 PRINT "No Integer solution from 1 to 100“ 
140 PRINT I; "Trials" : END 


Girls and boys problem 
Program is S-L-O-W. Be patient. 
40 Girls 50 Boys 
3950 Trials 


The computer program uses two FOR loops (Statements 40-120 and 50- 
1 10) to try every combination of values for B and G between 1 and 100 until a 
solution is found or until the program runs out of values. The variable I 
(Statement 60) is a counter which records the number of trials required to 
reach a solution. 
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The program is straightforward and finds a solution after 3950 trials. How- 
ever, it would have been a simple matter to substitute the value of B from the 
first equation in the second one and quickly solve the problem by hand or 
with the aid of a calculator. It is important to recognize that if a problem can 
easily be solved by other methods, the computer offers little or no advantage. 

Try this problem. You may or may not want to use your computer. If Mat- 
thew can beat Jeff by one-tenth of a mile in a two-mile race and Jeff can beat 
Steven by one-fifth of a mile in a two-mile race, by what distance could Mat- 
thew beat Steven in a two-mile race? (Hint: the answer is not 3/10 mile.) 


<*, 



33 


Brown's Books 


The use of a trial and error approach can generally be improved signifi- 
cantly if the combinations to be tried can be narrowed down in some way. 
The solution to this problem illustrates how the speed of obtaining a solution 
can be improved well over 100 fold by combining equations and eliminating 
certain solution possibilities. 

Here is the problem. Brown sold 48 books at a flea market, some for $3 
each, some for $5 each and others for $8. He collected a total of $175. He 
remembered having an even number of $5 books. Can you determine how 
many of each kind of book he had? 

The equations for solution are (letting T equal the number of $3 books, F 
the number of $5 books, and E the number of $8 books): 




3*T + 5*F + 8*E = 175 


The first program was written simply to try all possible combinations of T, 
F, and E from 1 to 48. It yields three solutions for the problem, although the 
two solutions with an odd number of $5 books can be eliminated leaving just 
the one desired solution. 

10 PRINT "Brown's Books" 

20 PRINT "L-O-N-G run time. Be patient." 

30 PRINTsPRINT "Computes all solutions." 

40 PRINT " $3 $5 $8* s PRINT 

50 FOR T=1 TO 48 

60 FOR F=1 TO 48 

70 FOR E=1 TO 48 

80 IF (T+F+E ) <>48 THEN 110 

90 IF T*3+F*5+E*8<> 175 THEN 110 

100 PRINT T;TAB(4)jF}TAB(10);E 

110 NEXT E 

120 NEXT F 

130 NEXT T 


8rown's Books 

L-O-N-G run tine. Be patient. 

Coiputes all solutions. 

$3 $5 $8 

34 13 1 

37 8 3 

40 3 5 


This program took approximately 22 minutes and 13 seconds to run on the 
Model 100 computer. A typical minicomputer (PDP-8/e) could run this 
problem in about 7.3 seconds. In either case, this is a long time to tie up the 
computer. 


34 


It is rather easy to combine the two equations into one by solving for T. 
The single equation is then: 

2*F + 5*E = 31 

In this equation, the limits can be reduced (from 48 used in the first run) 
since F cannot possibly be greater than 31/2 or 15.5 and E cannot be greater 
than 31/5 or 6.2. Making the appropriate program modifications leads to the 
second program. 

20 

50 

60 FOR F = 1 TO 15 

70 FOR E=1 TO 6 

80 IF ( F*2+E*5) <>31 THEN 110 

90 T=48-F-E 
130 

Brown's Books 

Computes all solutions. 

S3 $5 $8 

40 3 5 

37 8 3 

34 13 1 

Using this program produces a dramatic improvement in the time to solu- 
tion. On the Model 100, the time is approximately 2.5 seconds and on the 
PDP-8, about 0.16 seconds. 

Since the problem states that F must be even, a final modification which 
steps F by two in Statement 20, can be made. This version of the program 
takes only 1.25 seconds to run on the Model 100 and 0.06 seconds to run on 
the PDP-8. 


30 PRINT "Solutions for F even" 
60 FOR F=2 TO 14 STEP 2 


Brown's Books 
Solutions for F even 
$3 $5 $8 

37 8 3 


Notice the enormous improvement in computing time required for a solu- 
tion, over 1000 fold on the Model 100 and 100 fold on a PDP-8. Brute force 
certainly is inefficient! It is generally worthwhile to think through most prob- 
lems, particularly big ones, before rushing to the computer. The computer 
may be fast, but we just improved its performance by 1000 times by using a 
little common sense. 
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Intersection of Sets 


Two sets of numbers can be combined to yield a third set by the operation 
of intersection. The intersection of two sets A and B is the set that contains all 
elements that belong to both A and B. It does not contain any other elements. 
The intersection is usually written An B. 

For example if M = < 0,2,4, 6 > and K ^ < 1,2, 3, 4 >, then M D K = 
< 2,4 >. 

This program finds the intersection of two sets of numbers. It has been 
written to find the intersection of the two repetitive sets described in State- 
ments 30 and 40. In the sample run, Statement 30 describes the set 
x = < 1,3,5, ... 19 > and Statement 40 describes the set 
y = < 2,5,8, . . . 29 > . 

Notice that successive values of x increase by 2 and y by 3. 


10 PRINT "Sets X and Y Intersect at points:" 

30 FOR X=1 TO 19 STEP 2 

40 FOR Y = 2 TO 29 STEP 3 

50 IF X= Y THEN 90 

60 NEXT Y 

70 NEXT X 

80 END 

90 PRINT X; 

100 GOTO 70 


Sets X and Y Intersect at points: 
5 11 17 


However, if the set cannot be so neatly described, it may be desirable to 
rewrite the program to examine any set of data. This is done with the read 
statement which reads into x the data points in the data statement. The pro- 
gram is set up to use the same y set as the first program, but the x set is de- 
fined in the data statement. 

You should be able to see from the first combination of sets that if there is a 
numerical pattern in the sets which intersect, then there is also a pattern in 
the resulting intersecting set. In the example, the x values increase by 2 and 
the y values by 3, hence the values in the intersecting set increase by 2 x 3 = 
6. Although the intersection of these sets could easily be calculated by hand, 
the computer can be an aid in evaluating more complicated sets. 
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10 PRINT "Sets X and Y Intersect at points:" 
30 READ X 

40 FOR Y = 2 TO 29 STEP 3 

50 IF X=Y THEN 90 

60 NEXT Y 

70 GOTO 30 

80 ENO 

90 PRINT Xj 

100 GOTO 70 

200 DATA 2,3,8,9,14,15,20,21,26,27 


Sets X and Y Intersect at points: 
2 8 14 20 26 

?OD Error In 30 
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Prime Factors 

A prime factor is a positive integer that has no factor except itself and one. 
The first ten prime factors (or numbers) are 2, 3, 5, 7, 11, 13, 17, 19, 23, and 
29. The definition gives the basic method for determining whether a number 
is prime: divide by all smaller integers down to 2, testing whether the remain- 
der is zero for at least one of them. If not, the number is prime. 

But this is highly inefficient. It is obvious that a number is not prime if it is 
any even number greater than 2; hence only odd divisors need to be tried. 
Also, it is not necessary to try divisors greater than the square root of the 

number. 

Since the division method is inefficient, various schemes have been devised 
to avoid division. The basic idea underlying all such schemes is called the 
sieve of Eratosthenes (276 B.C.-195 B.C.). Imagine a list of odd numbers 
from 3 up. Strike out every third number after 3, every fifth number after 5, 
and so on. This will leave only prime numbers. 


10 CLS s PRINT "Prime factors. Enter 0 to stop. 
20 DIM A(100) 

30 C = 0 
40 PRINT 

50 INPUT "Your number" } M 
60 IF M=0 THEN END 
70 N=M 


80 X=0 

90 IF M>0 THEN 110 
100 PRINT "Invalid. 
110 1=1 
120 1 = 1+1 

130 IF DM THEN 190 
140 IF M/I<>INT(M/I) 


Again please. 


THEN 120 


150 X=X+1 

160 A(X)=I 

170 M=M/I 

180 GOTO 140 

190 IF X=1 THEN 250 

200 PRINT "Factors: * j Nj 

210 FOR L= 1 TO X 

220 PRINT A(L) ; 

230 NEXT L 

240 GOTO 260 

250 PRINT N" 1 s prime. "; 

260 PRINT : GOTO 40 


: GOTO 50 


H 
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Prlue factors. Enter 0 to stop. 

Your number? 105 
Factors: 105 357 

Your number? 72 

Factors: 72 2 2 2 3 3 

Your number? 365 
Factors: 365 5 73 

Your number? 89 
69 Is prime. 

Your number? 0 
Ok 


The program here finds the prime factors of any integer, or prints out “N is 
prime” if the integer has no proper divisors. 

Run this program for a large number of different integers and see if you 
can discover relationships between numbers and their prime factors. You 
should also try to figure out the method employed in the program to find the 
prime factors of any integer. To do this, you might want to draw a flowchart 
to show what is happening in the program. This will help you see the method 
used to find a prime factor and might help you in writing a program to gen- 
erate primes. 

In writing a program to generate prime factors, you can use the sieve 
method. However, as the numbers become very large, you will have to figure 
out a way to represent integers with more digits than your computer can han- 
dle at one time. (One approach is described on pp. 19-21 of Computers in 
Mathematics. ) 

Goldbach was a mathematician who made a conjecture that every even 
number greater than 4 can be written as the sum of two prime numbers (16 = 
11 + 5, 30 = 17 + 13, etc.). No one has ever proved it but no one has dis- 
proved it either. That is why it is called a conjecture. Can you write a pro- 
gram that will prove or disprove this conjecture? Or how about writing a 
program to prove Goldbach’s conjecture for even numbers up to 50? You 
should be able to write this program with 12 or fewer statements. 

Here is another problem involving prime numbers. Assume a life span of 
80 years. In what year of the 20th century (1900-1999) would a person have 
to be born to have the maximum number of birthdays occurring in prime 
years? The minimum number? 
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Greatest Common Divisor 


The greatest common divisor of a set of numbers is, as its name implies, the 
greatest integer that will divide into a set of two or more numbers. For exam- 
ple, the set of numbers 12, 20, and 28 have a greatest common divisor of 4. 
Nothing larger than 4 will divide evenly into all three numbers. 

This program will find the greatest common divisor for any set of integers. 
To run it, you simply input the number of integers in your set, type them in 
when requested and let the program calculate the GCD. The heart of the 
calculation is in Statement 150. 

Do you know the meaning of a relatively prime set of numbers? Can you 
figure out the meaning from the third sample run of the program or from 
runs of your own? How is a set of relatively prime numbers different from a 
set of prime factors? Can you find a set of 10 integers that is relatively prime? 


10 CLS-.PRINT "Greatest Common Divisor" 

20 DIM X{ 50 ) 

30 PRINT 

40 INPUT "Numbers In set" ; N 
50 PRINT "Enter numbers:" 

60 S=lE+25 
70 FOR K=1 TO N 
80 INPUT X(K) 

90 IF X( K ) >S THEN 110 
100 S=X ( K ) 

110 NEXT K 
120 G=0 

130 FOR M=2 TO 5 
140 FOR 1=1 TO N 

158 IF X(I)/M<>INT(X(I)/M) THEN 180 

160 NEXT I 

170 G=M 

180 NEXT M 

190 IF G> 0 THEN 220 

200 PRINT "Numbers are relatively prime. 

210 GOTO 230 

220 PRINT "GCD 1s" ; G 

230 INPUT "Another (Y or N)" ; Z$ 

240 IF Z$=“Y" OR Z$= "y " THEN 30 
250 PRINT "Okay, Bye for now." : END 


Greatest Common Divisor 

Numbers In set? 3 
Enter numbers: 

? 12 
? 36 
? 96 

GCD Is 4 

Another (Y or N)? y 

Numbers In set? 3 
Enter numbers: 

? 20 
? 36 
? 96 

GCD Is 4 

Another (Y or N ) ? y 

Numbers In set? 3 
Enter numbers: 

? 20 
? 36 
? 97 

Numbers are relatively prime 
Another (Y or N ) ? n 
Okay, Bye for now. 


In the last section we discussed prime numbers. Here is an interesting chal- 
lenge for you involving prime numbers. Until late 1982, the longest pro- 
gression of prime numbers in which all differed by the same number was 17. 
Prof. Paul Pritchard in the computer science department at Cornell Univer- 
sity wrote a program to determine if there was a longer progression. Using a 
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DEC VAX- 1 1/780, he found the string of 18 numbers shown below. He also 
discovered fourteen other 17-number progressions and ten 18-number pro- 
gressions, but none yet with 19 numbers. He believes there is at least one; can 
you find it? 


107928278317 

117851061187 

127773844057 

137696626927 

147619409797 

157542192667 

167464975537 

177387758407 

187310541277 


197233324147 

207156107017 

217078889887 

227001672757 

236924455627 

246847238497 

256770021367 

266692804237 

276615587107 
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Cryptarithmetic Problems 

Cryptarithmetic or alphametic problems are arithmetic expressions in 
which the digits are replaced by letters of the alphabet. Each digit is asso- 
ciated with a letter to produce an interesting statement, for example: 

SEND 
+ MORE 

MONEY 

If the college student who sent this message to his father needed $106.52 
for plane fare home, this was the right message to send since this combination 
of letters has one unique solution, in particular: 

9567 
-f 1085 

10652 

However, if the student let things go to the last moment and was in more of 
a rush, he might have reworded the message: 

WIRE 
+ MOR E 
MONEY 

In this case, how much should his dad send? Earlier in the book, trial and 
error approaches to solving problems were discussed. It was noted that the 
brute force approach of trying every alternative was sometimes appropriate. 
Is it in this case? 

No! The number of possible alternative solutions is the factorial of the 
number of different letters in the alphametic expression, i.e., 8! or 40,320. A 
program to try out every one of these possibilities would run for a 1-o-n-g 
time. 

In this case it is much more efficient to apply some common sense to nar- 
row down the number of alternatives. The best approach to this process is to 
divide up the search space into large classes (or sets), according to a common 
property shared by members of each class, and then attempt to eliminate en- 
tire classes by the method of contradiction. 

Consider the “WIRE + MORE = MONEY” problem. Can E = 0? Since 
E + E = Y, Y must also equal 0, contradicting the fact that Y and E must be 
different digits. Thus, the entire class of solutions in which E = 0 can be 
ruled out. 

Consider E = 3. Now Y = 6 and there is no carry to the next column. So 
in this column R + R — E or E + 10 if a carry is involved. But in either case 
E must be an even number since 2R is always even; this contradicts the 
assumption that E = 3. 

By following this type of classificatory contradiction process for each of the 
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digits in the order E, R, I, O and N, the computer program will search out all 
possible solutions to the problem. Unlike the "send more money” problem, 
the ’’wire more money” problem has five possible solutions. A smart father 
would choose the lowest solution and wire his son $103.48. 


10 CLS 

20 PRINT "Cryptarl thmetl c Solver" 

30 PRINT " WIRE" 

40 PRINT " + M 0 R E" 

50 PRINT ■ " 

60 PRINT "MONEY" 

70 M=1 

60 FOR E=2 TO 9 
90 Y=E+E 

100 IF Y>=10 THEN 130 

110 C1 = 0 

120 GOTO 150 

130 Cl = l 

140 Y=Y-10 

150 FOR R = 0 TO 9 

160 IF R=M OR R=E OR R=Y THEN 450 

170 IF R+R+C1=E THEN 200 

180 IF R+R+C1=E+10 THEN 220 

190 GOTO 450 

200 C2=0 

210 GOTO 230 

220 C2 = l 

230 FOR 1=1 TO 9 

240 IF I=M OR I=E OR I=Y OR I=R THEN 440 
250 FOR 0=0 TO 9 

260 IF 0=M OR 0=E OR 0=Y OR 0=R OR 0=1 THEN 430 

270 N=I+0+C2 

280 IF N> = 10 THEN 310 

290 C3=0 

300 GOTO 330 

310 C3=l 

320 N=N-10 

330 IF N=M OR N=E OR N=Y OR N=R OR N=I OR N=0 THEN 430 
340 FOR W=0 TO 9 

350 IF W=M OR W=E OR W=Y OR W=R OR W=I OR W=0 OR W=N THEN 420 
360 IF O+10OW+M+C3 THEN 420 
370 PRINT 

380 PRINT " " ; W; I ; R ; E 

390 PRINT " + " ; M ; 0 ; R ; E 

400 PRINT " " 

410 PRINT M; 0; N; E j Y 

420 NEXT W 

430 NEXT 0 

440 NEXT I 

450 NEXT R 

460 NEXT E 

470 END 
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Cryptarl thuetl c Solver 



W 

I 

R 

E 

+ 

M 

0 

R 

E 

M 

0 

N 

E 

Y 


9 

7 

6 

2 

+ 

1 

0 

6 

2 

1 

0 

8 

2 

4 


9 

2 

7 

4 

+ 

1 

0 

7 

4 

1 

0 

3 

4 

8 



9 

5 

7 

4 

+ 

1 

0 

7 

4 

1 

0 

6 

4 

8 


9 

2 

8 

7 

+ 

1 

0 

8 

7 

1 

0 

3 

7 

4 


9 

5 

8 

7 

+ 

1 

0 

8 

7 

1 

0 

6 

7 

4 


There are other approaches to solving cryptarithmetic problems, but all of 
them benefit greatly from reducing the search space as much as possible be- 
fore putting the problem on the computer. See if you can devise another 
successful approach and write a program to implement it. 

Here are some problems for you to try. 


DONALD 
+ GERALD 

TWO 
x TWO 

ABCDE 
x 4 

ROBERT 

THREE 

EDCBA 

ABC 

THE 

SPRING 

x DE 

EARTH 

RAINS 

FFC 

VENUS 

BRING 

1 I— 'V-' 

DEC 

HGBC 

SATURN 
4- URANUS 

NEPTUNE 

4- GREEN 
PLAINS 

ONE 

FORTY 

FIVE 

TWO 

TEN 

-FOUR 

4- FIVE 

4- TEN 

ONE 

EIGHT 

SIXTY 

+ ONE 
TWO 

VIOLIN + VIOLIN + VIOLA 4- CELLO 

= QUARTET 

THREE + NINE = 

EIGHT 4- FOUR 
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Sailors and Monkey Problem 


There are many variations of the sailors and monkey problem. Here is one 
of them. 

Five sailors and a monkey were on an island. One evening the sailors 
rounded up all the coconuts they could find and put them in a large pile. Be- 
ing exhausted from working so hard, they decided to wait and divide them up 
equally in the morning. During the night, a sailor awoke and separated the 
nuts into five equal piles, but had one nut left over which he gave to the mon- 
key. He took one pile, hid it, and pushed the other four together and went 
back to sleep. He was followed in this action by the other four sailors, each of 
whom did exactly the same thing. Next morning the remaining nuts were di- 
vided equally with one remaining nut going to the monkey. What is the small- 
est number of coconuts with which they could have begun? 

Although there is an elegant algebraic solution to this problem, a more 
suitable approach for the computer is that of working backwards. A typical 
solution to a problem can be thought of as a path that leads from the given 
information to the goal. However, in this case the goal, or final state, is 
known, thus it is easier to start there and work backwards to the initial state. 

As mentioned at the outset, many sailor and monkey problems exist, in 
fact, an infinite number of them. For example, instead of five sailors, there 
could be three or six or 14. Thus it is desirable to devise a general solution in- 
stead of just one to solve one specific problem. 

In the flowchart and computer program, S is the number of sailors and A is 
the number of coconuts that each sailor received in the final division of the 
pile. Since one coconut was given to the monkey at each division, the total 
number of coconuts left in the morning must be S x A + 1. But this pile came 
from pushing together S - 1 equal piles. Thus, the key condition that must 
hold for (S x A + 1) / (S - 1) to be an integer K, which represents the num- 
ber of coconuts that the last sailor stole from a pile of S x K + 1 coconuts. 
But this pile is the result of pushing together S - 1 equal piles by the previous 
thief, so again (S x K + 1) / (S - 1) is an integer and so on back through all S 
raids on the pile. 

Note in the flowchart (and program) that the first trial value for A is 1 
(Statement 40). In Statement 80 this value is increased by 1 until the value of 
(S x K + 1)/ (S - 1) is an integer as tested for in Statement 70. This process is 
then continued until the counter for the second loop, N (nighttime pile di- 
visions) equals the number of sailors. 

Although the program will work for any number of sailors, it takes a fairly 
long time to run for more than five. Remembering what you have learned ear- 
lier in this chapter, can you devise a way to make the program more efficient? 
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10 CIS : 
30 INPUT 
40 A=1 
50 K=A 
60 N=0 


PRINT "Sailors and Monkey." 
"Number of sa11ors" } S 


H J F < S * k+ 1 ) / (S-1)=INT((S*K+1)/(S-1)) THEN 100 
80 A=A+l 

90 GOTO 50 
100 N=N+l 

110 K= ( S*K+1 ) / ( S-l ) 

120 IF N=S THEN 140 
136 GOTO 70 

140 PRINT "Total coconuts =";S*K+1 

}50 " In the "owning, each sailor gets" 

160 PRINT A; "coconuts . " 


Sailors and Monkey. 

Nuiber of sailors? 3 
Total coconuts = 79 
In the lornlng, each sailor gets 
7 coconuts. 

Sailors and Monkey. 

Number of sailors? 5 
Total coconuts = 15621 
In the i»orn1nc» ( each sailor gets 
1023 coconuts. 
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Super Accuracy 

Under normal circumstances, your computer performs computations to six 
or seven digits of accuracy. Double precision computations increase accuracy 
to 13 digits or so. 

However, it is possible to do computations one digit at a time and assign 
each digit to an element in an array. This will achieve virtually any desired 
accuracy. Any, up to the maximum array size that is. 

This example program performs the rather simple operation of successively 
doubling a number (which is the same as raising 2 to a power). 

If the number to be represented is 8192, then: 

A(4) A(3) A(2) A(l) 

8 1 * 9 2 


To add this to itself, first the rightmost digits are added: A(l) + A(l). 

If there is a carry, the variable C is set equal to 1, otherwise C is 0. The total 
is put into B in Line 100. 

If B is less than 10, there is no carry (C=0) and the new A(l) equals B. If 
B is greater than 10 there is a carry (C= 1) and the new A(l) equals B-10. 

This operation is continued for all the digits (D) of the number and, when 
it is finished, the new number is printed in Lines 190-210. 

If A(N) is printed followed by a semicolon(;) for tight packing, the Basic 
print routine would leave a space in front of each digit (for the sign) and a 
space after each digit (for readability). In the program here, these spaces are 
not wanted, hence the print routine in Line 200 is used which prints the string 
value of the ASCII value of each digit (which is the same as the digit itself) 
but without the spaces. 

This program, incidentally solves the challenge to calculate the number of 
moves in the Towers of Brahma problem (see “Change for Any Amount to 
$5.00). The approach is also used in the next section, “Palindromes.” 


10 PRINT "Computes 2 to the Nth to any ac 
26 OIM A(100) 

30 M=0 s C=0 
40 0=0 

50 FOR 1 = 1 TO 100 : A ( I ) = 0 : NEXT I 
60 A ( 1 ) = 1 
70 1=0 : C=0 
80 M=M+1 
90 1=1+1 

100 B=A( I )+A( I )+C 

110 IF B<10 THEN C=8 : GOTO 140 


120 B=B-10 

130 C = 1 

140 A ( I ) =8 

150 IF KD THEN 90 

160 IF C=1 THEN 90 

170 PRINT M; 


180 

190 

200 

210 

220 

230 


curacy" s PRINT 


D=I 

FOR N= I TO 1 STEP -1 
PRINT CHRS(A(N)+48)i 
NEXT N 
PRINT 
GOTO 70 
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Computes 2 to the Nth to 

any accuracy 

1 

2 

67 

147573952589676412928 

2 

4 

68 

295147905179352825856 

3 

8 

69 

590295810358705651712 

4 

16 

70 

1180591620717411303424 

5 

32 

71 

2361183241434822606848 

6 

64 

72 

4722366482869645213696 

7 

128 

73 

9444732965739290427392 

8 

256 

74 

18889465931478580854784 

9 

512 

75 

37778931862957161709568 

10 

1024 

76 

75557863725914323419136 

11 

2048 

77 

151115727451828646838272 

12 

4096 

78 

302231454903657293676544 

13 

8192 

79 

604462909807314587353088 

14 

16384 

80 

1208925819614629174706176 

15 

32768 

81 

2417851639229258349412352 

16 

65536 

82 

4835703278458516698824704 

17 

131072 

83 

9671406556917033397649408 

18 

262144 

84 

19342813113834066795298816 

19 

524288 

85 

38685626227668133590597632 

20 

1048576 

86 

773712524*55336267181195264 

21 

2097152 

87 

154742504910672534362390528 

22 

4194304 

88 

309485009821345068724781056 

23 

8388608 

89 

618970019642690137449562112 

24 

16777216 

90 

1237940039285380274899124224 

25 

33554432 

91 

2475880078570760549798248448 

26 

67108864 

92 

4951760157141521099596496896 

27 

134217728 

93 

9903520314283042199192993792 

28 

268435456 

94 

19807040628566084398385987584 

29 

536870912 

95 

39614081257132168796771975168 

30 

1073741824 

96 

79228162514264337593543950336 

31 

2147483648 

97 

158456325028528675187087900672 

32 

4294967296 

98 

316912650057057350374175801344 

33 

8589934592 

99 

633825300114114700748351602688 

34 

17179869184 

100 

1267650600228229401496703205376 

35 

34359738368 

101 

2535301200456458802993406410752 

36 

68719476736 

102 

5070602400912917605986812821504 

37 

137438953472 

103 

10141204801825835211973625643008 

38 

274877906944 

104 

20282409603651670423947251286016 

39 

549755813888 

105 

40564819207303340847894502572032 

40 

1099511627776 

106 

81129638414606681695789005144064 

41 

2199023255552 

107 

162259276829213363391578010288128 

42 

4398046511104 

108 

324518553658426726783156020576256 

43 

8796093022208 

109 

649037107316853453566312041152512 

44 

17592186044416 

110 

1298074214633706907132624082305024 

45 

35184372088832 

111 

2596148429267413814265248164610848 

46 

70368744177664 

112 

5192296858534827628530496329220096 

47 

140737488355328 

113 

10384593717069655257060992658440192 

48 

281474976710656 

114 

20769187434139310514121985316880384 

49 

562949953421312 

115 

41538374868278621628243970633760768 

50 

1125899906842624 

116 

83076749736557242056487941267521536 

51 

2251799813685248 

117 

166153499473114484112975882535043072 

52 

4503599627370496 

118 

332306998946228968225951765070086144 

53 

9007199254740992 

119 

664613997892457936451903530140172288 

54 

18014398509481984 

120 

1329227995784915872903807060280344576 

55 

36028797018963968 

121 

2658455991569831745807614120560689152 

56 

72057594037927936 

122 

5316911983139663491615228241121378304 

57 

144115188075855872 

123 

10633823966279326983230456482242756608 

58 

288230376151711744 

124 

21267647932558653966460912964485513216 

59 

576460752303423488 

125 42535295865117387932921825928971026432 

60 

1152921504606846976 

126 

85070591730234615865843651857942052864 

61 

2305843009213693952 

127 

170141183460469231731687303715884105728 

62 

4611686018427387904 

128 

340282366920938463463374607431768211456 

63 

9223372036854775808 

129 

680564733841876926926749214863536422912 

64 

18446744073709551616 

130 

1361129467683753853853498429727072845824 

65 

36893488147419103232 

131 

2722258935367507707706996859454145691648 

66 

73786976294838206464 

132 

5444517870735015415413993718908291383296 


Palindromes 


A palindrome is a word, verse, or number that reads the same backwards 
or forwards. For example, the words “mom” and “eye” are palindromes. So 
are each of the lines in this verse: 

Egad, a base life defiles a bad age 

Doom an evil deed, liven a mood 

Harass sensuousness, Sarah 

Golf; No, sir, prefer prison-flog 

Ban campus motto, “Bottoms up, MacNab” 

Numeric palindromes are those numbers which read the same backward as 
forward. The examination of these numbers is a field rich with possibilities 
for creative computing. 

One conjecture concerning palindromes raises an interesting unanswered 
question. Begin with any positive integer. If it is not a palindrome, reverse its 
digits and add the two numbers. If the sum is not a palindrome, treat it as the 
original number and continue. The process stops when a palindrome is ob- 
tained. For example, beginning with 78: 

78 
+ 87 

165 
+ 561 

726 
+ 627 

1353 
+ 3531 

4884 

The conjecture, often assumed true, is that this process will always lead to 
a palindrome. And indeed that is just what usually happens. Most numbers 
less than 10,000 will produce a palindrome in less than 24 additions. But 
there is a real thorn in the side of this conjecture, the number 196. Can you 
determine if a palindrome will ever be produced with a starting number of 
196? 

The number 196 will produce 1675 after two reversals, but after 100 rever- 
sals the resultant sum has 47 digits and is still not palindromic. Why mention 
1675? Because ten other numbers under 1000 will also lead to the sum of 
1675 and thus may not become palindromic. The first five of these numbers 
are 196, 295, 394, 493, and 592. What are the other five? 
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The program here will accept any number as a starting value and complete 
the process of adding the successive reversals and testing if the sum is a palin- 
drome. Try it with some numbers and see if you can identify any patterns. 


10 PRINT "Tests conjecture that adding reversed" 
PRINT "digits will produce a palindrome." 


30 DIM B ( 50 ) 

40 PRINT 

50 INPUT "Start number";A 
60 E=0 
70 E=E+1 
60 A=A/ 10 

90 IF INT ( A ) > 0 THEN 70 
100 FOR C = E TO 1 STEP -1 
110 A=A*10 

120 B(C) = INT(A-10*INT(A/10) ) 

130 NEXT C 
140 0=0 

150 FOR C = 1 TO I NT ( E / 2 ) 

160 IF B( C ) =B( E+l-C ) THEN 180 
170 D=1 
180 NEXT C 

190 FOR C=E TO 1 STEP -1 

200 PRINT CHR$(B(C)+48)i 

210 NEXT C 

220 IF D=1 THEN 260 

230 PRINT " Palindrome! "j 

240 PRINT 

250 GOTO 40 

260 PRINT " Not yet" 

270 IF E/2<=INT (E/2) THEN 290 
280 B(INT(E/2)+l)=2*B(INT(E/2)+l) 
290 FOR C=1 TO INT( E/2 ) 

300 B(C)=B(C)+B(E+1-C) 

310 NEXT C 

320 FOR C=1 TO INT( E/2 ) 

330 B( E+l-C ) =B{ C ) 

340 NEXT C 

350 B ( E+l ) = 0 

360 FOR C=1 TO E 

370 B(C+1)=B(C+1)+INT(B(C)/10) 

380 B(C)=B(C)-10*INT(B(C)/10) 

390 NEXT C 

400 IF B( E+l ) < = 0 THEN 140 
410 E=E+1 
420 GOTO 140 


Tests conjecture that adding reversed 
digits will produce a palindrome. 

Start number? 19 
19 Not yet 
110 Not yet 
121 Palindrome! 

Start number? 35 
35 Not yet 
88 Palindrome! 


Tests conjecture that adding reversed 
digits will produce a palindrome. 

Start number? 89 
89 Not yet 
187 Not yet 
968 Not yet 
1837 Not yet 
9218 Not yet 
17347 Not yet 
91718 Not yet 
173437 Not yet 
907808 Not yet 
1716517 Not yet 
8872688 Not yet 
17735476 Not yet 
85189247 Not yet 
159487405 Not yet 
664272356 Not yet 
1317544822 Not yet 
3602001953 Not yet 
7193004016 Not yet 
13297007933 Not yet 
47267887164 Not yet 
93445163436 Not yet 
176881317877 Not yet 
955594586548 Not yet 
1801200002187 Not yet 
8813200023188 Palindrome! 

Start number? 196 
196 Not yet 
887 Not yet 
1675 Not yet 
7436 Not yet 
13783 Not yet 
52514 Not yet 
94039 Not yet 
187088 Not yet 
1067869 Not yet 
18755470 Not yet 
18211171 Not yet 
35322452 Not yet 
68744805 Not yet 
111589511 Not yet 
227574622 Not yet 
454050344 Not yet 
897100798 Not yet 
1794102596 Not yet 
8746117567 Not yet 
16403234045 Not yet 
70446464506 Not yet 
130992928913 Not yet 
450822227944 Not yet 
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960544455998 Not yet 
1800098901607 Not yet 
8801197801088 Not yet 
17602285712176 Not yet 
84724043932847 Not yet 
159547977975595 Not yet 
755127757721546 Not yet 
1400255515443103 Not yet 
4413700670963144 Not yet 
8827391431036288 Not yet 
17653692772973576 Not yet 
85191620502609247 Not yet 
159482241005228405 Not yet 
664364741147513356 Not yet 
1317626482294916822 Not yet 
3603815405135183953 Not yet 
7197630720180367016 Not yet 
13305261530450734933 Not yet 


47248966933966985264 Not yet 
93507933867933969538 Not yet 
177104867844767940077 Not yet 
947154635293536341848 Not yet 
1795298270686072793597 Not yet 
9749270977546801719568 Not yet 
18408442064004592449047 Not yet 
92502871604050616929528 Not yet 
175095833209091234750057 Not yet 
925153265399993573340628 Not yet 
1751196640799987135692157 Not yet 
9264161958699957602603728 Not yet 
17537224026299926194218357 Not yet 
92918473189299188236491928 Not yet 
175837936477498486373973857 Not yet 
934217310162393261013712428 Not yet 
1758434620324786522027424867 Not yet 
9442681822581660752291773438 Not yet 


Using this method, write a program that examines all the integers between 
1 and 10,000 excluding those that sum to 1675 at any point. What does this 
show? By the way, you will have to devise a way to deal with 14-digit integers 
which are larger than your computer can normally handle. 

Huh? Is this program “too hot to hoot?” 
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Convergence and Recursion 


The computer is especially suitable for doing repetitive and tedious calcula- 
tions. Two mathematical approaches for solving problems that involve repet- 
itive calculations are convergence and recursion. 

Some problems can be reasonably easily stated in words or described with a 
few simple equations but there are many possible solutions. For example, how 
many ways can you make change for a dime? It is simply stated and the num- 
ber of ways can be enumerated fairly easily: two nickels, one nickel and five 
pennies, or ten pennies, three ways in all. But if you want to solve for all the 
ways of making change for a dollar or five dollars, it would be nice to have 
some help. 

Help on this kind of problem comes from a class of computer program that 
simply breaks the problem into smaller ones and counts up all the alternative 
solutions according to a set of rules. But an even more powerful technique is 
known as recursion. Using this technique, a simple solving algorithm or rou- 
tine is set up to solve the smallest subset of the problem. The unique power in 
a recursive routine comes from the ability of the routine actuallly being able 
to call itself. This is discussed further in the second program in this section. 

Another approach for solving problems that do not have an exact answer is 
that of successive approximations. For example, the exact value of pi, e or the 
length of an irregular curve cannot be precisely determined. But by means of 
increasingly accurate approximations, it is possible to approach the desired 
value from above or below or to converge on it from two directions. The last 
four programs in this chapter illustrate successive approximations and 
convergence. 


Change For a Dollar 

Even though there is very little you can buy for a penny these days, the 
coin will probably be around for some time to come since it is needed to make 
change for odd amounts of sales tax and to fill up penny collections. 

Today U.S. coinage consists of five coins: penny, nickel, dime, quarter, and 
half dollar. How many ways can coins of these denominations be used to 
make change for one dollar? For example, one way is two half dollars, an- 
other is one half dollar and two quarters, and so on. Make a best guess now 
and write it down before you read further. 

There are several different ways to approach a problem of this kind. One is 
to break it down into smaller, more easily solved problems. In other words, 
how many ways can you make change for a quarter? For a dime? You would 
solve these subproblems and combine the answers to give the overall solution. 

If you were more mathematically inclined, you could write a series of equa- 
tions relating each piece of change to every other one and to the dollar and 
solve them. 

A third approach is to do the problem by writing down combinations until 
all the different possibilities are exhausted (or until you are exhausted) and 
then count them all up. This might be called solving the problem by exhaus- 
tion and is a method quite suitable for putting on the computer. 

Write a program that uses this approach to solve the problem. If you use 
loops and count by one, it could take a long time for the computer to run 
through all the possible combinations, possibly many hours. 

Also, if you want to print out all the possible combinations, be warned that 
the printing could also take quite some time and a fair amount of paper. 
There are more combinations than you might think! 

In fact, most people will not be able to guess the answer to this problem, or 
even come close. Ask several of your friends how many ways they think a 
dollar can be changed. Record all the responses and then tabulate them on 
your computer. What is the mean (average) of all the guesses? The extremes? 

The program included here uses the first method discussed to solve the 
problem, in particular, breaking down the problem into subproblems and 
then combining the solutions into one final answer. 

First, the main problem is broken into the next smaller one of making 
change for half dollars. There are three such problems: no half dollars (H = 
0), one half dollar (H = 1), and two half dollars (H = 2). The last problem is 
trivial since there is only one way, but the other two need to be broken down 
further. 

This is done by dividing the remaining money into quarters and consid- 
ering the subproblems on down to the lower denominations. As the number 
of subproblems is expanded, each one becomes easier to solve. In fact, 
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subgoals, which can be solved in only one way, are finally reached. For exam- 
ple, if H = 1, Q = 1, D = 2, and N = 0, then the pennies (P) must equal 
five in order that the total equal 100. 

Notice that at the quarter, dime and nickel stages, adjustments are made in 
the limits of the loops depending upon how much money there is left to 
change. For example, if H = 1, the only possible subgoals for quarters are 0, 
1, and 2, but not 3 or 4. Also notice that there is no need to test combinations 
of coins to see if they add up to 100, nor is it necessary to include the penny 
as a variable. Simply counting the number of subgoals is sufficient since each 
one can be solved in only one way. 


10 C = 0 

20 FOR H = 0 TO 2 
30 FOR Q=0 TO (4-2*H) 

40 FOR 0=0 TO ( 10-5*H-2. 5*Q) 

50 FOR N = 0 TO (20-10*H-5*Q-2*D) 

60 C=C+1 
70 NEXT N 
80 NEXT D 
90 NEXT Q 
100 NEXT H 

110 PRINT "A dollar can be changed In* 
120 PRINT Cj "different ways." 


A dollar can be changed In 
292 different ways. 


Try to make some changes in this program or write a new one to solve the 
following problems. Say you want two quarters in your change to play some 
video games. In how many ways can a dollar be changed to provide at least 
two quarters? 

Visiting a small town, you find the parking meters still take pennies. In 
how many ways might you get change so that you had at least three pennies? 
Is this any different than the number of ways that would give you five pen- 
nies? Say you want to make a phone call also; in how many ways can you 
change a dollar to produce at least four pennies and one dime? 
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Change For Any Amount to $5.00 

Another way to attack the change problem in the previous section is by 
means of the programming technique called recursion. Get familiar with this 
one — it is very powerful! Donald Piele and Larry Wood described this 
method in an issue of Creative Computing. 

First, define the variables which represent the number of ways to make 
change for n cents using the coins specified: 

A Only pennies 
B Nickels and pennies 
C Dimes, nickels, and pennies 
D Quarters, dimes, nickels, and pennies 
E Halves, quarters, dimes, nickels, and pennies. 

Initially, there are two subproblems in making change for n cents. In the 
first, no half dollars are used, and D is the number of ways to change n cents. 
Second, when one or more halves are used, after one is paid, there remain 
n-50 cents to pay which can be done in E n _ 50 ways. 

Since these two cases are mutually exclusive, it can be inferred that E n = 
D n + En_ 5 o. Similarly, 

D n = C n + D n _25 
Cn = B n + C n _io 

B n = A n + B n _s 

Now, begin with the simplest case and build up to Ejoo- First of all, it is 
easy to understand why E 0 = 1. From above, when n = 50, E go = D 50 + 
E 0 , and it is possible to make change for 50 cents only one more way if half 
dollars are allowed. Therefore E 0 = 1. Likewise, D 0 = C Q = B 0 =A 0 = 1. 
It is also true that A„ = 1 for all values of n since there is only one way to 
make change using only pennies. Now the recursive relationships can be used 
to solve the original problem. 

This is the strategy used in the program. It also has the added advantage 
that it can count the number of ways of making change (with coins) for any 
specified amount. 

Now it is your turn. Can you modify the program here to include dollar 
bills so it could count the number of ways to make change for any amount up 
to $10.00? 

Using any method of change making you prefer, write a program to make 
change for one ruble. Russian coins come in denominations of 1, 2, 3, 5, 10, 
15, 20, 50, and 100 kopecks. There are 100 kopecks in one ruble. 
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10 DIM A(ldl) t B(iei) l C(ltl) l D(iei),E(iei) 

20 PRINT: INPUT "Amount to be changed"; X 

30 IF X>5 THEN PRINT "Too much; again please..." ; GOTO 20 
40 M=INT (20*X)+1 

50 A(l) = ls B(l)=lsC(l)=l: D(l) = l: E(l)=l 
60 FOR 0=2 TO M 
70 A( 0) =1 

80 B(0)=A(0)+B( J-l ) 

90 C( J)=B(0) 

100 IF 0<=2 THEN 120 
110 C(0)=B(0)+C(0-2) 

120 D(J)=C(J) 

130 IF J<=5 THEN 150 
140 0( J)=C( J)+D( J-5) 

150 E(J)=D(J) 

160 IF O<=10 THEN 180 
170 E( J)*D( J)+E( J-10) 

180 NEXT J 

190 PRINT "You can make change for $";X 
200 PRINT "1n";E(M) { "different ways." 

210 GOTO 20 


Amount to be changed? .1 
You can make change for $ .1 
In 4 different ways. 

Amount to be changed? .12 
You can make change for $ .12 
In 4 different ways. 

Amount to be changed? 1 
You can make change for $ 1 
In 292 different ways. 

Amount to be changed? 3.89 
You can make change for $ 3.89 
In 23124 different ways. 

Amount to be changed? 5.19 
Too much; again please... 

Amount to be changed? 4.19 
You can make change for $ 4.19 
In 30270 different ways. 
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Perhaps the most famous problem used to demonstrate the principles of 
recursion is the Towers of Brahma. It is sometimes called the Towers of Ha- 
noi or Pharoah’s Needles. Here is the problem in as close to original form as 
possible. You should be able to solve it with a relatively short program using 
recursion. 

In the great temple at Benares beneath the dome which marks the center of 
the world rests a brass plate in which are fixed three diamond needles, each a 
cubit high and as thick as the body of a bee. On one of these needles, at the 
Creation, God placed 64 discs of pure gold, the largest disc resting on the 
brass plate and the others getting smaller and smaller up to the top one. This 
is the Tower of Brahma. 

Day and night unceasingly, the priests transfer the discs from one needle to 
another, according to the fixed and immutable laws of Brahma. These laws 
require that the priest on duty must not move more than one disc at a time 
and that he must place this disc on a needle so there is no smaller disc below 
it. When the 64 discs shall have been thus transferred from the needle which, 
at the Creation, God placed them, to one of the other needles; tower, temple, 
and Brahmans alike will crumble into dust, and with a thunderclap, the 
world will vanish. 

If the priests were to effect one transfer every second, and work 24 hours 
per day for each day of the year, it would take them 58,454,204,609 decades 
plus slightly more than six years to perform the feat, assuming they never 
made a mistake — for one small slip would undo all the work. 

How many transfers are required to fulfill the prophecy? Try out your pro- 
gram with fewer discs than 64 to make sure you are on the right track. Here 
is a table of the first few transfers: 


Discs 

Moves 

Discs 

Moves 

1 

1 

6 

63 

2 

3 

7 

127 

3 

7 

8 

255 

4 

15 

9 

511 

5 

31 

10 

1023 
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Converge on e and Pi 


An incredibly important mathematical constant is designated by the small 
letter e. This constant is both irrational and transcendental. Look up those 
terms in a dictionary or math book if you wish, or just plunge on to the next 
paragraph. 

The constant e was first derived by John Napier, also the inventor of loga- 
rithms, to whom we owe an eternal debt of gratitude. Why? If e had never 
been discovered, advances in mathematics, physics, and astronomy would 
have lagged a century or more, because e is the base of all natural logarithms 
and these logarithms are the basis for many branches of science and 
mathematics. 

How is e calculated? The constant e is the limiting value of this expression 
as n approaches infinity: 

e = (l + 1 A) n 


Its exact value can never be found, but to 15 places e equals 
2.718281828459045 . . . How can e be calculated? First take 1 -% and square 
it, that equals 2- 1 /,. Then cube 1 - l / 3 and you get 2.3686. Raising l- l / 4 to the 
fourth power gives 2.414, and so on. Write a program for this method and 
have it print out the initial value of e and each value after each additional 
fraction is added on. 

Another approach is to expand the expression above using the binominal 
theorem and, again, letting n approach infinity. The expression for the expan- 
sion is: 


e = 1 + 1 


^ 2! T 3! T 4! ^ 5! " ‘ 


n\ 


Now here is where the computer can again be of some assistance. Since 3! is 
3*2! and 4! is 4*3!, all the calculations need not be done for each additional 
fraction. Look at the program and particularly note the calculations in State- 
ments 70 and 90. 


10 

PRINT "Converge on e" 

Converge 

20 

PRINT “with binominal series." 

2 

30 

E=1 : 1=0 

2.5 

40 

I-I + l 

2.66667 

50 

0=1 

2.70833 

60 

FOR J=1 TO I 

2.71667 

70 

0=0*0 

2.71806 

80 

NEXT 0 

2.71825 

90 

E=E+1/D 

2.71828 

100 

PRINT E 

2.71828 

110 

GOTO 40 

2.71828 

2.71828 
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Pi is another important mathematical constant that is irrational (meaning 
its exact value can never be determined) and transcendental (meaning it is not 
the solution to any algebraic equation). Interestingly, pi was known and used 
by the ancients. Archimedes, who lived in the second century B.C., by using a 
regular polygon of 96 sides (nearly a circle), proved that the value of pi was 
less than 22/7 and greater than 3 - 10 / 71 , a remarkable achievement for the 
mathematics of his day. 

Ptolemy in 150 A.D. used the value of 3.1416 for pi and in the middle of 
the sixteenth century the amazing fraction 355/113 was discovered, giving 
the value of pi accurately to six decimal places. 

Incidentally, in 1897, the General Assembly of Indiana passed a bill ruling 

that the value of pi was four. 

Several infinite series can be used to grind out increasingly accurate values 
for pi. One such series is (l-Va + % - Vr + %•■■)■ This series is called an 
arithmetic series and converges very slowly. The program here displays only 
every 500th value of the series. Don’t be alarmed if the program does not 
seem to be running very fast; a fair amount of calculating is going on between 
each value that is printed. 


Converge on pi 
using arithmetic series 
3.1436 
3. 14059 
3.14226 
3. 14109 
3. 142 
3. 14126 
3.14188 
3. 14135 
3.14182 

3. 1414 
3.14178 
3.14144 
3.14176 
3.14146 
3. 14174 

3. 14148 
3. 14172 

3. 14149 
3. 14171 

3. 1415 
3. 1417 


10 PRINT "Converge on pi” 

20 PRINT "using arithmetic series" 
30 S=1 : 1=1 = Q=8 : P=0 
40 Q=Q+1 
50 P=P+S/I 
60 1=1+2 
70 S=-S 

80 IF Q<499 THEN 40 
90 Q=0 

100 PRINT P*4 
110 GOTO 40 


Actually, the approach used by Archimedes converges much more quickly 
and, with the aid of a computer, it is possible to go far beyond the 96-sided 
polygon used by Archimedes. 
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B C 



His approach was to construct inscribed and circumscribed polygons and 
measure the perimeters to approximate the circumference of a circle. Con- 
sider a polygon circumscribed around a circle of radius 1. The perimeter 
equals the length of one side times the number of sides. Since the tangent of x 
= AB/BC, but BC = 1, then tan (x) = AB and the length of a side = 2 tan 
(x). Since the circumference = 2 pi r and r = 1, then pi is the circumference 
(or perimeter of an n-sided polygon) divided by 2. 

Similar trigonometry leads to the perimeter of an inscribed polygon being 
equal to the number of sides times sin (x) * cos (x). 

The second program produces values for pi using inscribed and circum- 
scribed polygons. Unfortunately, there is one large flaw in the program, be- 
cause degrees must be converted into radians in Statement 60. This means, of 
course, that you must already know the value of pi, since the conversion fac- 
tor is 360 degrees divided by 2 pi. 

Setting this flaw aside, it is interesting to note how quickly this program 
converges on the value of pi compared to the preceding one. That is because 
this one converges geometrically rather than arithmetically. 

Can you figure out a geometric convergence to the value of pi that does not 
require that you know it (or a conversion factor) before you start? 


10 PRINT "Converge on pi using polygons." 

20 PRINT : PRINT "Sides Inscribed Cf rcumscrl bed" 


30 C=360/57. 2958 
40 N=6 
50 N = 2*N 
60 X=C/N 
70 PRINT N; 

80 PRINT TAB(7);N*SIN(X)*C0S(X)/2; 
90 PRINT TAB(20) ; N*TAN(X)/2 
100 GOTO 50 


Converge on pi using polygons. 


SI des 

Inscrl bed 

Circumscribed 

12 

2.59808 

3.4641 

24 

3 

3.21539 

48 

3. 10583 

3. 15966 

96 

3. 13263 

3. 14609 

192 

3. 13935 

3. 14271 

384 

3. 14103 

3. 14187 

768 

3. 14145 

3.14166 

1536 

3. 14156 

3. 14161 

3872 

3. 14158 

3. 1416 

6144 

3.14159 

3. 14159 

12288 

3. 14159 

3. 14159 
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Convergence on Pi Revisited 

In answer to the question posed in the last paragraph of the preceding sec- 
tion, here is a way to converge on pi without knowing its value beforehand. 

As in the previous program, the basic approach is to add up the length of 
the sides on an inscribed polygon and divide by 2r to obtain a value for pi. 
The program starts with a square (four sides) and doubles the number of sides 
each time. If the old side length is S, then the length of S’ of a side of a new 
polygon with twice as many sides is obtained by applying the Pythagorean 
theorem. In particular, 

X 2 + (S/2) 2 = R 2 
(R-X) 2 + (S/2) 2 = (S’) 2 

Thus, — — — — — — ^ 

S’ = V (R-n/R 2 - (S/2) 2 ) 2 + (S/2) 2 

It is easy to reduce this formula algebraically, but accuracy suffers if this is 
done. Also, you will find that S*S is slightly more accurate than ST 2. 

Unfortunately, striving for maximum accuracy is somewhat moot on a 
computer that does not have double precision arithmetic. Notice that ac- 
curacy does not improve with more than 1024 sides and, indeed, as numbers 
in the calculations start to exceed the capacity of the computer (4194304 
sides), the accuracy starts to deteriorate badly. 

There is yet another method to compute pi by convergence. It uses discov- 
eries of Gregory and Euler. Gregory discovered the formula for arctangent: 
Arctan x = x - x 3 /3 + x 5 /5 - x 7 /7 + ... 

Euler came up with a rather interesting formula for pi: 
pi = 4 (arctan (V 2 ) + arctan (Vs)) 

See if you can combine these two formulas to calculate pi. If you are very 
clever, you can do your calculation to yield far more than the seven decimal 
place accuracy obtained by the programs presented here so far. 



of Circle 
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10 PRINT “Converge on pi with* 

20 PRINT "Inscribed polygons." 

30 PRINT : PRINT "Sides Perimeter" 
40 R=10 
50 Z = 2 
60 N=Z+Z 
70 S=R*SQR(Z) 

80 FOR K=1 TO 28 

90 PRINT Nj TAB( 10 ) ; (S*N)/(Z*R) 

100 Y=S*S/ (Z+Z) 

110 X=R-SQR(R*R-Y) 

120 S=SQR(X*X+Y) 

130 N=N*Z 
140 NEXT K 


Converge on pi with 
Inscribed polygons. 


SI des 

Perl meter 

4 

2.82843 

8 

3.06147 

16 

3.12145 

32 

3.13655 

64 

3.14033 

128 

3. 14128 

256 

3.14151 

512 

3.14157 

1024 

3. 14159 

2048 

3.14159 

4096 

3. 14159 

8192 

3.14159 

16384 

3. 14159 

32768 

3.14159 

65536 

3. 14159 

131072 

3. 14159 

262144 

3.14159 

524288 

3.14159 
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Length of Any Curve 

The previous programs have demonstrated how it is possible to compute a 
very accurate value of pi by adding together the length of the sides of a poly- 
gon as it approaches a circle and dividing by 2r. 

Using a similar approach, it should be possible to inscribe a polygon, or 
portions of a polygon, inside any regular curve and thus determine the length 
of the curve. This program approximates the length of any curve as defined in 
Statement 260 by dividing it into an increasing number of subintervals and 
computing the sum of the secants (a straight line that cuts a curve at two or 
more points). 

To run the program, you must enter the formula or equation describing 
your curve in Statement 260 in the form: 

260 FN = your function of X 

You then type RUN and enter the end points of the curve you want to use 
in your calculation. These points are entered in the form of the abscissa, 
which means the horizontal (or y) coordinate of the point. 

The program is written to sum the successive secant lengths and to cal- 
culate the percent of change in each summation compared to the preceding 
one. The sample run uses the function 2x 3 + 3x 2 -2x + 3. Note the substan- 
tial improvements in the length calculation as the number of intervals in- 
creases from 2 to 16, but the rather slight improvements beyond 16. 

Try this program with different curves and functions. It might help to plot 
the function first (remember the program to do that?) and then compute its 
length. Is this method of length calculation more accurate for a function with 
ntion with one or more changes? Why? 
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10 SCREEN 0,0 

20 PRINT "Length of any curve. " 

30 PRINT "Define function in Subroutine 260." 
40 INPUT "Min. & max. abscissa (y)";P,Q 
50 PRINT "Intvls Secant Length X Change" 

60 S1=0 

70 FOR N= 1 TO 9 
80 E=2~ ( N-l ) 

90 H= { Q-P ) / E 
100 S=0 

110 FOR 1=0 TO E-l 
120 X=P+I*H+H: GOSUB 260 
130 F1=FN 

140 X=P+I*H: GOSUB 260 
150 F2=FN 

160 L=SQR((F1-F2)"2+H*H) 

170 S=S+L 
180 NEXT I 

190 IF S1>0 THEN 220 

200 PRINT E;TAB(5)jS;TAB(20) ; "No prev value" 
210 GOTO 240 

220 PC=((ABS(S1-S))/S1)*100 
230 PRINT EjTAB(5)jS;TAB(20);PC 
240 S1=S 
250 NEXT N 

260 FN = 2*X / '3+3*X A 2-2*X+3 
270 RETURN 


Length of any cur 
Define function 1 
Min. & max. abscl 
Intvls Secant Le 


n Subroutine 260. 
ssa ( y ) ? - 1 , 6 
ngth X Change 


1 

525.047 

No prev value 

2 

525. 158 

.0212849 

4 

529.652 

.855722 

8 

531.017 

.257703 

16 

531.964 

. 178376 

32 

532.017 

9.82135E-03 

64 

532. 042 

4. 73811E-03 

128 

532.049 

1 . 28485E-03 


260 FN=SIN( X) 


Length of any cur 
Define function 1 
Min. & max. abscl 
Intvls Secant Le 


ve. 

n Subroutine 260. 
ssa ( y ) ? - 1 , 1 
ngth X Change 


1 

2.61387 

No prev value 

2 

2.61387 

9. 1213E-06 

4 

2.62005 

.236643 

8 

2.62214 

.0797049 

16 

2.6227 

.0212401 

32 

2. 62284 

5. 4089E-03 

64 

2.62287 

1. 35442E-03 

128 

2.62288 

3. 09059E-04 


Converge on a Square Root 


Since most square roots are irrational, methods used to calculate them usu- 
ally involve successive approximations. Although you can simply call up the 
square root function in Basic or on many pocket calculators, it is interesting 
to explore various methods of calculating square roots without these built-in 
functions. After all, these built-in functions are nothing more than successive 
approximation routines already installed in the machine. 

Obviously, a square root is the inverse of the operation of squaring a num- 
ber. All of the methods of calculating square roots use this fact, but the way 
in which it is employed is quite different in various calculators and 
computers. 

The program here calculates an upper and lower limit for the square root 
of a number and, by successive approximations, pinches the root to within a 
smaller and smaller interval until it reaches the desired level of accuracy. 

The starting value for the lower limit is 0 and for the upper limit the num- 
ber, Z, whose square root is sought. The program then divides this interval 
into ten steps by simply dividing the difference between the numbers by 10. 
The variable I is increased from the lower limit to the upper one by the value 
of the step, S. At any point, if I squared becomes greater than Z, a new upper 
limit is set to I and a new lower limit is set to I - S. 

This method converges very quickly and adds approximately one decimal 
place of accuracy with each pass beyond the third. What happens when you 
enter into the program a number that has an exact square root such as 25 or 
49? Why? 

Another approach to calculating square roots by successive approxima- 
tions is to start with a trial root, X. If X * X is less than the original number 

N, then increase the trial value by a 0.1. If X * X is greater than N, return to 
the previous value. This is the first digit of the root. Now, start advancing by 

O. 01. Continuing in this way, one digit is developed at a time until the desired 
precision is reached. 

This method is quite suitable and fast for square roots of numbers less than 
1. A good first trial root value is 0.1. But is it suitable for larger numbers? 
How should a starting trial value be determined? In this method, especially 
for numbers greater than 10, the initial trial value for the root matters a great 
deal in determining the length of time it will take for the calculation to con- 
verge. Write a program using this method and compare the speed and ac- 
curacy with the program in the book. 
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10 PRINT "Square Root Convergence" 
20 PRINT : INPUT "Your number";Z 
30 E = 1 E - 0 5 

40 PRINT : PRINT "Lower Upper" 
50 PRINT "Limit Limit" 

60 A=0 : B=Z 

70 S= ( B-A ) / 10 

80 PRINT A; TAB( 18 ) ; 6 

90 IF ABS(A*B-Z)<E THEN 180 

100 FOR I=A TO B STEP S 

110 IF Z<I*I THEN 156 

120 NEXT I 

130 B=B*1 0 

140 GOTO 70 

150 B= I 

160 A=I-S 

170 GOTO 70 

180 PRINT "Accurate to .00001" 

190 PRINT "Average =" ; (A+B)/2 


Square Root Convergence Square Root Convergence 

Your number? 3 Your number? 54 


Lower 

Upper 

Lower 

Upper 

Limit 

Limit 

Limit 

Limit 

0 

3 

0 

54 

1.5 

1.8 

5.4 

10. 8 

1.71 

1.74 

7.02 

7.56 

1.731 

1.734 

7. 344 

7.398 

1.7319 

1.7322 

7. 344 

7.3494 

Accurate 

to .00001 

7.34832 

7.34886 

Average = 

» 1.73205 

7.34843 

7.34848 



7.34847 

7.34847 



Accurate 

to .00001 

Square Root Convergence 

Average « 

= 7.34847 

Your number? 49 





Square Root Convergence 

Lower 

Upper 



Limit 

Limit 

Your number? 25 

0 

49 



4.9 

9.8 

Lower 

Upper 

6.86 

7.35 

Limit 

Limit 

6.958 

7.007 

0 

25 

6.9972 

7.0021 

5 

7.5 

6.99965 

7.00014 

5 

5.25 

7 

7.00005 

5 

5.025 

7 

7 

5 

5.0025 

7 

7 

5 

5.00025 

Accurate 

to .00001 

5 

5.00003 

Average = 

7 

5 

5 



5 

5 



Accurate 

to .00001 



Average = 

5 
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Compounding 


As with successive approximations and recursion, compounding requires 
many repetitive calculations. Some compound interest and growth situations 
can be represented by a formula, but for many problems, solving by repetitive 
calculations is an excellent approach. 

There is nothing magical about compounding; you generally start with an 
initial amount of money, number of animals, etc. This quantity then grows or 
diminishes by a certain percentage at set intervals. The new amount is the old 
amount increased or decreased by this percentage. Repeat this calculation 
over and over again, and you have solved the problem. 

Five different types of problems which involve compounding of some sort 
are described in this chapter. Try the extra problems that follow some of the 
programs; you may be surprised at the results. 
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Indians and Interest 


Here is a simple compound interest problem that produces an astonishing 
answer. The problem has to do with the sale of Manhattan Island to the 
Dutch for approximately $24 worth of trinkets and beads. 

If the $24 that the Indians received in 1626 had been deposited in a bank 
paying 5- 3 / 4 % interest compounded annually, how much would it amount to 
in 1983? 

The program here solves this little problem by making use of the formula 
to calculate compound interest. In particular, if P dollars are invested at an 
interest rate of R (expressed as a decimal) and compounded N times, then the 
total amount A is given by the formula: 

A = P(1 + R) N 

How much was gained in 1983 alone? How much was gained in the decade 
from 1974 to 1983? You can change the value of N in Statement 30 to get the 
answers to these questions. But a better way might be to enter the ending year 
with an input statement. 

Can you read a number expressed in the E (exponential format)? In a more 
conventional format, the number is $11,176,500,000 or $11.2 billion. 

There are many other ways you can improve the program and make it 
more suitable for general purpose compound interest calculations. Modify it 
to accept any starting and ending year, any rate of interest, and any starting 
principal amount. 

The problem as stated is somewhat unrealistic since banks were not paying 
interest rates of 5- 3 / 4 % from 1626 to 1983. Change the program to calculate 
the total amount based on the following interest rates: 


Years 

1626-1830 

1831-1870 

1871-1910 

1911-1921 

1922-1929 

1930-1940 

1941-1945 

1946-1960 

1961-1980 

1981-1983 


Interest Rate 
1.5% 
2 . 0 % 
3.0% 
3.5% 
6.5% 
2.3% 
3.5% 
5.3% 
6.5% 
9.5% 


70 


10 P = 24 
20 R=. 0575 
30 N=1983-1626 
40 A=P*( ( 1+R ) *N ) 

50 PRINT "Indians got $24 for Manhattan In" 
60 PRINT "1626. At 5.75%, by 1983 It would" 
70 PRINT "be worth $"iA 


Indians got $24 for Manhattan In 
1626. At 5.75%, by 1983 It would 
be worth $ 1.11765E+10 
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Systematic Savings 


In the previous program, compound interest was calculated by means of a 
formula well known to bankers, money lenders and real estate agents. How- 
ever, if you did not know the formula, how would you calculate interest on 
money in a savings account by hand or with a calculator? 

You would probably begin by multiplying the principal amount P by the 
interest rate R and adding that amount to the original principal at the start of 
the second year. Doing that for all the years the money is in the bank will 
yield a final amount. Why not write a program to perform the calculations in 
this manner rather than use a formula? Here is such a program. 

It is set up for an initial principal amount of 100 (Statement 50), an interest 
rate of 10% (R = 0.1 in Statement 60) and ten years (N = 10 in Statement 
70). Naturally these values could be read in using input statements. The 
clever calculation in Statement 100 rounds off the amount to two decimal 
places (dollars and cents). 

In contrast with the program in the previous section, this one does not use 
a compound interest formula, but simply adds the interest each year to the 
growing principal amount in Statement 90. 


10 CLS : 

PRINT "Calculates Interest on $100" : 
PRINT "Invested at 10% for 10 years." : 


PRINT 

30 PRINT "At end Current" 

40 PRINT "of year balance" 

50 P = 1 00 
60 R = . 1 
70 N=10 

80 FOR 1=1 TO N 

90 p=p+p*R 

100 PRINT I;TAB(10); INT(P*100+.5)/100 
110 NEXT I 


Calculates Interest 

i nvested 

at 10% for 

At end 

Current 

of year 

bal ance 

1 

110 

2 

121 

3 

133. 1 

4 

146.41 

5 

161.05 

6 

177.16 

7 

194.87 

8 

214.36 

9 

235.79 

10 

259.37 


on 

10 


$100 

years. 


Now, how could this program be modified to allow for a plan of systematic 
saving? In other words, instead of letting the $100 lie around all lonely while 
it is compounding, each year you add another $100 to it. With this program, 
making the modification for systematic savings is easy: Statement 105 is 
added to add the new deposit each year to the ever growing principal. 
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10 CLS : 

PRINT "Calculates Interest on $100" : 

PRINT "Invested at 10% every year for 10 yrs " • 
PRINT 

30 PRINT "End Current" 

40 PRINT "year Invested Balance" 

50 P= 1 00 
55 C=100 
60 R = . 1 
70 N= 10 

80 FOR 1=1 TO N 
90 P=P+P*R 

100 PRINT I;TAB(6);I*C;TAB(15);INT(P*100+.5)/100 
105 P=P+C 
110 NEXT I 


Calculates Interest on $100 
Invested at 10% every year for 10 yrs. 

End Current 

year Invested Balance 


1 

100 

110 

2 

200 

231 

3 

300 

364. 1 

4 

400 

510.51 

5 

500 

671.56 

6 

600 

848.72 

7 

700 

1043. 59 

8 

800 

1257.95 

9 

900 

1493.74 

10 

1000 

1753. 12 


In an effort to attract depositors, many banks over the last 20 years have 
started offering savings and investment accounts in which the deposits are 
compounded more often than annually. In the early 50’s, many banks went to 
quarterly compounding; in the late 50’s, to daily; and in the early 60’s, some 
“competitive” S&L’s went to continuous compounding. 

How much does more frequent compounding really mean? Modify either 
of the two programs to compute the interest for more frequent compounding. 
What is the difference in interest for one year for annual, quarterly, monthly, 
daily, and continuous (every second) compounding on a principal amount of 
$1000 invested at 8%? How about for ten years? 

Incidentally, if you want to use the formula method in the previous section 
for this calculation, the formula for P principal invested at R rate com- 
pounded N times per year is: 

A = P(1 + R/N) N 

Try this problem which uses the same principles of compounding. Con- 
sumer prices rose an average of 8.8% in 1980. While the government keeps 
trying to bring inflation under control, they don’t seem to be meeting with 
much success. Assuming that prices continue to go up this much (8.8%) ev- 
ery year, how much will a $6000 economy car (1980 dollars) cost in the year 
2000? How much will it cost when you are 65 years old? 
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Systematic Savings Revisited 

Using the formula for compound interest and systematic savings, this pro- 
gram will calculate the amount accumulated after a given period of time. 

’ When you run the program, it will ask how much you wish to save each 
month, the number of compounding periods in a year, the interest rate, and 
the length of time you wish to continue your systematic savings program. The 
program will then calculate the total amount at the end of that period. 

From the preceding programs, you should be able to see that a systematic 
savings program is a very effective way to accumulate a nest egg for the 
future. 


10 PRINT "Calculates Interest and balance for" : 
PRINT "a systematic savings program." : 
PRINT 

20 INPUT "Monthly depos1t";A 
30 INPUT "Compounding periods/year" B 
40 INPUT "Interest rate (xx.x)";R 
50 R=R/100 

60 INPUT "Number of years'jN 
70 PRINT ; 

PRINT "End" f TAB(14); "Current" 

80 PRINT "Year Invested Balance" 

90 C=12*A 
100 P*12*A 
110 FOR 1=1 TO N 
120 P=P*( 1+R/B)~B 

130 PRINT I}TAB(6);I*C{TAB(14)jINT(P*100+.5)/100 
140 P=P+C 
150 NEXT I 


Calculates Interest and balance for 
a systematic savings program. 

Monthly deposit? 10 
Compounding periods/year? 1 
Interest rate (xx.x)? 8.5 
Number of years? 8 

End Current 

Year Invested Balance 

1 120 130.2 

2 240 271.47 

3 360 424.74 

4 480 591.04 

5 600 771.48 

6 720 967.26 

7 840 1179.68 

8 960 1410.15 


Calculates Interest and balance for 
a systematic savings program. 

Monthly deposit? 10 
Compounding periods/year? 365 
Interest rate (xx.x)? 8.5 
Number of years? 8 


End Current 


Year 

Invested 

Balance 

1 

120 

130.64 

2 

240 

272.88 

3 

360 

427.73 

4 

480 

596.32 

5 

600 

779.86 

6 

720 

979.68 

7 

840 

1197. 23 

8 

960 

1434.08 


Calculates Interest and balance for 
a systematic savings program. 

Monthly deposit? 200 
Compounding perl ods / year? 12 
Interest rate (xx.x)? 12 
Number of years? 5 

End Current 

Year Invested Balance 


1 

2400 

2704. 38 

2 

4800 

5751.74 

3 

7200 

9185.57 

4 

9600 

13054. 9 

5 

12000 

17415 


Try to combine the things that you have learned from the programs in 
these sections to write a completely generalized program for systematic sav- 
ings. It should accept the following information as input: 

• Initial deposit 

• Frequency of periodic deposits 

• Amount of each deposit 

• Interest rate 

• Interest compounding frequency 

• Length of time of savings program 

Your program should produce output in tabular form showing years (or 
other time periods), amount invested without interest, total amount with in- 
terest, and the interest alone. 
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Loan Payments 


Although saving money in a systematic way is a noble goal, many people 
frequently find themselves talking to a different bank officer, namely the one 
in charge of loans. 

This program will calculate the payments for a loan for a period of one 
year or longer. The program asks you to enter the key facets of the loan in 
question: amount borrowed, annual rate of interest, interval between pay- 
ments and term of the loan in years. 

The sample run shows a relatively short-term loan (2 years) of $3000 at a 
bargain basement interest rate of 8.5%: The monthly payment of the loan is 
computed to be $136.37 and the total interest $272.79. Why is the total in- 
terest not $255 (8.5% of $3000)? Instead it is 9.09% Is this a mistake in the 
program? 

Try the program for some longer term loans at real world interest rates. 
For example, run it for an automobile loan of $8000 at 12% over a 5-year 
period. Perhaps you can see from this that systematically saving your money 
and paying cash for an item makes more sense than time payments. 

Run the program to calculate the mortgage payments on a $150,000 house 
with a $40,000 down payment. Try it with a 16% interest rate stretched over 
a 30-year period. Ouch! Look at all that interest! 

At the heart of this program are Statements 110-140. Can you see what is 
happening in these calculations? 


The second version of the program uses a PRINT USING statement in 
Line 50 to produce somewhat nicer looking output. 


10 PRINT "Calculates payaent schedule on a loan." : 
PRINT 

20 INPUT *A»ount borrowed";A 
30 INPUT "Interest rate (xx.x)";I 
40 INPUT "Payaent Interval (ionths)*;P 
50 INPUT "Ter* (years)" } Y 

60 INPUT "Output tables (1) or totals (2) only'jV 
70 PRINT 

80 IF V= 2 THEN 110 

90 PRINT "Table Is In *jP;"aonth Intervals" : 

PRINT 

100 PRINT "Per- Interest Principal" : 

PRINT "lod Due Due" s 

PRINT 

110 Z=(Y*12)/P 

120 K=( I*(P/12) ) / 1 00 

130 E»A*K/(1-1/(1+KPZ) 

140 E=INT(E*100+.5)/100 
150 C=A : 

F=0 s 
D1 = 0 : 

T1 = 0 
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160 T1=T1+1 

170 IF T1>Z THEN 250 

180 B=T1 : 

C=C-F s 
D=C*K : 

F = E-D 

190 D=INT(D*100+. 5)/100 

200 F=INT(F*100+. 5)/100 

210 D1=D1+D 

220 IF V = 2 THEN 160 

230 PRINT B;TAB(5);D;TAB(14)jF 

240 GOTO 160 

250 IF V=1 THEN 300 

260 D1 = INT(D1*100+. 5)/100 

270 PRINT "Principal $"iA 

280 PRINT "Interest " 5 D1 

290 GOTO 320 

300 PRINT 

310 PRINT "Ttl " ; D1 ; TAB( 14 ) ; A 
320 X=INT((Dl+A)*100+.5)/100 
330 X1=X/((Y*12)/P) 

340 X1=INT(100*X1+. 5)/100 
358 PRINT "Total " ; X 

360 PRINT "Each payment $"iXl 


Calculates payment schedule on a loan. 


Amount borrowed? 3000 
Interest rate (xx.x)? 8.5 
Payment Interval (months)? 1 
Term (years)? 2 

Output tables (1) or totals (2) 


only? 


2 


Principal $ 3000 
Interest 272.79 

Total 3272.79 

Each payment $ 136.37 


Calculates payment schedule on a loan. 


Amount borrowed? 8000 
Interest rate (xx.x)? 13 
Payment Interval (months)? 3 
Term (years)? 5 

Output tables (1) or totals (2) 


only? 


1 


Table Is In 3 month Intervals 


Per- 

Interest 

Prl ncl pal 

iod 

Due 

Due 

1 

260 

290. 23 

2 

250.57 

299.66 

3 

240.83 

309.4 

4 

230.77 

319.46 

5 

220.39 

329.84 

6 

209.67 

340.56 

7 

198.6 

351.63 

8 

187.17 

363. 06 

9 

175.38 

374.85 

10 

163.19 

387.04 

11 

150.61 

399.62 

12 

137.63 

412.6 

13 

124.22 

426.01 

14 

110.37 

439.86 

15 

96.08 

454. 15 

16 

81.32 

468.91 

17 

66. 08 

484. 15 

18 

50.34 

499. 89 

19 

34. 1 

516. 13 

20 

17.32 

532.91 

Ttl 

3004.64 

8000 

Tota 

1 

11004.6 

Each 

payment 

$ 558.23 
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Interest on Credit Purchases 

All too frequently, the rate of interest to be charged on a loan or credit pur- 
chase is hidden in the small type. After all, the bank or car dealer or finance 
company wants to convince you that you can afford that new car or home 
improvement or whatever. 

This program calculates the interest rate on a loan given the principal value 
of the loan, the number of payments, and the amount per payment. To make 
things even easier, the program will accept the cash purchase price of the arti- 
cle and the down payment and automatically compute the principal value of 
the loan. 

The sample run shows a rather unrealistic loan for an item costing $88.99. 
The down payment was $10.00 and the loan is over a period of 18 months; 
each monthly payment is $4.85. The program run shows that the actual in- 
terest rate is a modest 5.69%. 

Now try the program with a more realistic loan. A termite company in 
New Jersey advertises a total treatment for your house for only $200; just 
$29.95 down and 24 monthly payments of $11.95 per month. Is this the bar- 
gain that it seems? What is the annual rate of interest? 
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10 PRINT "Calculates Interest rate of* 

PRINT " a purchase on cred 1t." s 

20 INPUT "Purchase pr1ce":P 
30 INPUT "Down payment":D 
40 INPUT "Number of payments":N 
50 INPUT "Payments per month";M 
60 INPUT "Amount per payment" ; A 
76 P1=P-D 
80 T=A*N-P1 
90 Y*N/(M*12) 

100 R=INT(10000*T/P1/Y+. 51/100 
110 PRINT : 

PRINT "Interest rate ts'jRj"** 


Calculates Interest rate of 
a purchase on credit. 

Purchase price? 88.99 
Down payment? 10 
Nuiber of payments? 18 
Payments per month? 1 
Amount per payment? 4.85 

Interest rate Is 7.01 » 


Calculates Interest rate of 
a purchase on credit. 

Purchase price? 1000 
Down payment? 0 
Number of payments? 24 
Payments per month? 2 
Amount per payment? 50 

Interest rate Is 20 % 


Purchase price? 1000 
Down payment? 0 
Number of payments? 24 
Payments per month? 1 
Amount per payment? 50 

Interest rate Is 10 I 
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Population Growth 


So far, all the compounding problems in this chapter have involved 
money — interest, savings, and loans. But many fascinating compounding 
problems do not involve money. Consider the following problem. 

In 1960, the population figures for the United States and Mexico were 180 
million and 85 million respectively. The annual population growth rate for 
the United States was 1.23% and for Mexico 2.23%. If these growth rates 
remain steady, in some distant year the population of Mexico will exceed that 
of the United States. In what year will this occur? 

The program uses the method of accumulating the principal amount rather 
than a formula, and applies it to both populations. The sample run reveals 
that the population of Mexico will overtake that of the U.S. in the not too 
distant future. 


ie 11=180000000# 

20 M=85000000# 

30 RU=1 . 0123 
40 RM=1 . 0223 
50 Y=1960 
60 U=U*RU 
70 M=M*RM 
80 Y=Y+1 

90 IF M<U THEN 60 

100 PRINT "Year = " 5 Y 

110 PRINT TAB(9) 5 "Population" 

120 PRINT "U.S. A. "jU 

130 PRINT "Mexico " : M 


Year = 2037 

Popul atl on 
U.S. A. 4 . 61403E+08 
Mexico 4. 64465E+08 


But what if the Mexican people were diligently trying to bring their 
increasing population under control and their annual growth rate was 
increasing by only 0.001% per year, compared to the increase of the U.S. 
growth rate of 0.01% per year. If that were the case, would the population of 
Mexico ever exceed that of the U.S.? 

Insert Statements 60, 70, and 85 into the first program, run it and find out 
the answer to the question. 


60 U=U*(1. 01*RU) 

70 M=M*( 1 . 001*RM ) 

85 IF Y > 9999 THEN 100 
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Going back to the first program, run it with the population data from the 
U.S. (180 million) and California (15.7 million) from 1960. At that time, the 
annual growth rates were 1.23% and 3.7% respectively. Running the pro- 
gram will indicate the year that the population of California will exceed that 
of the United States. This, of course, is nonsense. Where is the discrepancy? 

This program might make more sense if it were restated to ask in what year 
the population of California would exceed that of the remainder of the U.S., if 
ever? 

Compounding can also be used to solve other kinds of population growth 
problems. Consider the bristleworm. The bristleworm can reproduce by split- 
ting itself into 24 segments, each of which grows a new head and tail. What is 
the maximum number of bristleworms that could be obtained in this fashion, 
starting with only one worm, after ten splittings? Assuming a splitting occurs 
every 22 days, how many offspring will one worm produce in a year? When 
will the earth be overrun with bristleworms? 

Here is another problem for which the principle of compounding is useful. 
It takes nature about 500 years to produce one inch of topsoil. Many years 
ago, the United States had an average depth of almost nine inches of this 
good dirt, but as of 1975, the country was down to about six inches. This type 
of soil is necessary, of course, for growing food. 

Careless management of our soil causes about 1 % per year to erode away; 
it is then lost forever. Once soil depth reaches three inches or less, it is impos- 
sible to grow crops on a large scale. Write a program to calculate the year in 
which the U.S. will have less the 3" of topsoil, assuming that it continues to 
erode away at 1% per year. 

Will you be alive then? Will your children be alive? Will anyone be alive? 
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Probability 


Statistics and probability are subjects that bring up all kinds of images. 
Some people who have not had a pleasant time in a required college statistics 
course keep as far away from the subject as possible. To other people, statis- 
tics is something with which devious manufacturers can mislead you about 
their products. 

Not long ago, one foreign automobile manufacturer boasted that 90% of 
their cars sold in the United States in the last seven years were still on the 
road. This sounds like excellent reliability. But consider the fact that this 
manufacturer was rapidly expanding in the U.S. market and 65% of the cars 
they had sold in the U.S. had been sold in just the previous two years. You 
would expect that all of these would still be running. 

If the manufacturer had sold 10% of their 7-year volume in Years 1 to 3, 
and most of these cars were not running, then the advertising claim loses 
much of its meaning. 

When approached in a logical, step-by-step manner, statistics and probabil- 
ity are not difficult. In fact, they can be a great deal of fun. 

Some of the programs use a formula while others simulate the event in 
which we are interested. The results are the same, but you may find that the 
simulations help you to understand exactly what is happening. 
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Pascal's Triangle— Calculated 

Pascal’s Triangle is quite fascinating. As you can see from the portion of 
one reproduced below, each row is symmetrical. Each row also contains the 
coefficients for a binomial expansion. The sums across the ascending diag- 
onals form the Fibonacci sequence. The sums across the rows are all powers 
of two. Each row corresponds to the digits of a power of 1 1. Every element is 
the sum of the two above it. And, in case you care, all the elements in it are 
identities in combinatorial theory. 

The ways this marvelous triangle can be generated are as varied and in- 
teresting as its properties, though perhaps more difficult to figure out. Here is 
one way to use a computer to generate the triangle. 

Any element can be found be adding together the two elements immedi- 
ately above it. The program uses this principle to produce a triangle eight lev- 
els deep. Lines 10-30 set the rightmost diagonal to 1. Each element is stored 
in the two-dimensional variable P(R,C) with R denoting the row and C 
denoting the “column.” The variable T (Line 50) simply leaves some blank 
spaces so the output resembles a triangle. The crux of the calculation is in 
Line 70 in which the value of each new element is calculated. 

There is another interesting way to calculate Pascal’s Triangle in even 
fewer statements than the program here. It generates the triangle one element 
at a time and does not use any arrays or two-dimensional variables. Can you 
figure out how to write such a program? 


5 SCREEN 0,0 

10 FOR C= 1 TO 8 

20 P ( C , C ) =1 

30 NEXT C 

40 FOR R=1 TO 8 

50 T=20-R*2 

60 FOR C=2 TO R+l 

70 P(R+1,C)=P(R,C)+P{R,C-1) 

80 IF P(R, C)=0 THEN 110 

90 PRINT TAB(T)jP(R,C); 

100 T=T+4 
110 NEXT C 
120 PRINT 
130 NEXT R 
140 GOTO 140 


1 

1 1 
1 2 1 

13 3 1 

1 4 6 4 1 

1 5 10 10 5 1 

1 6 15 20 15 6 1 
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Pascal's Triangle — by Probability 

This program simulates the dropping of balls through a triangular array 
shown below. At each level, a ball is equally likely to fall either to the left or 
right. At the bottom of the array a cup is placed at each end point; these cups 
collect the balls. After each group of balls has been dropped, the number of 
balls in each cup is tallied and displayed. 


10 PRINT "Pascal's Triangle by Probability" 
20 INPUT “Balls to droD":M 
30 INPUT "Number of 1 
40 PRINT "Position 
50 FOR N= 1 TO M 
60 T = 0 

70 FOR 1=1 TO K 
80 IF RND( 1 )>. 5 THEN 
90 T=T+1 
100 NEXT L 
110 B( T+l ) =B(T+1 ) +1 
120 NEXT N 
130 FOR L=1 TO K+l 
140 PRINT L, B(L) 

150 B(L)=0 
160 NEXT L 


u u u u u u 




84 


Pascal ’ s Triangle by 
Balls to drop? 648 
Nimber of levels? 6 
Position frails 

1 11 

2 68 

3 146 

4 202 

5 150 

6 55 

7 8 


Probabl 1 1 ty 


The number of balls landing in the various cups at each level, when re- 
duced to the lowest common denominator, should approximate the numbers 
in Pascal’s Triangle. Do you know why they should? 

Try a few runs of the program with a different number of balls. Do the 
numbers obtained really approximate those in Pascal’s Triangle? If, at each 
dividing point in the array, every other ball went in the opposite direction, 
then Pascal’s Triangle would be produced. Try running the program with 4 
balls at Level 2; with 8 balls at Level 3; and with 16 balls at Level 4. Do your 
results look like those of the previous program which calculated the triangle? 
Does the approximation come closer as the number of balls is increased? 

How can you determine how close your results from this program are to 
the exact value of Pascal’s Triangle? One way is to take the number of balls 
that dropped into each cup on a given level and divide that by the total num- 
ber of balls divided by the theoretical sum of the row. So, in the sample run at 
Level 4, you would divide the balls in each cup by 62.5 (1000/16). You then 
compare that to the “correct” number and compute the percent difference. 

Here is the result of this procedure for Level 4 (actually the fifth row) in 
the sample run: 


Cup 

Balls 

+ 62.5 

Correct 

Deviation 

1 

52 

.83 

1 

17.00% 

2 

253 

4.05 

4 

1.25% 

3 

370 

5.92 

6 

1.33% 

4 

266 

4.26 

4 

6.50% 

5 

59 

.94 

1 

6.00% 


Why are the outside cups “off’ by a greater percentage than those closer to 
the center? Does this also happen on other levels with different numbers of 
balls? 
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Common Birthdays 


Here is an interesting little problem in probability. In a group of ten people 
selected at random, what is the probability that any of them will share the 
same birthday? How about a group of 20 people? Of 50 people? 

Conversely, how many people would you need in a group such that there is 
a 50% probability that at least two of them have the same birthday? How 
many people would be needed for a 90% probability of an overlap? 

Try to answer these questions, either by guess or by calculation, before you 
look at the output from the program. 

This program provides a painless introduction to the world of statistics. 
The calculation is actually quite trivial. The probability that any person in a 
group has a birthday on January 1 is 1/365. If our group has only two people 
in it, the probability that both of them have a birthday on January 1 is 1/365 
times 1/365, or a very small number indeed. The probability that they have a 
common birthday is 365 times the very small number just obtained, or about 
0.27%. 

However, if there are three people in the group, the probability goes up 
slightly. Call the three people Betsy, Ken, and Larry. From the reasoning 
above, we know that there is a 0.27% probability that Betsy and Ken have 
the same birthday; also a 0.27 % probability that Betsy and Larry share a 
birthday; and finally, a 0.27% chance that Larry and Ken have a common 
birthday. Hence, the total probability for a group of three people is three 
times the probability for just two people. 

A group of four increases the probability over that of two people by a fac- 
tor of six, five people by a factor of 10, six people by 15, seven people by 21, 
and so on. There is a progression here, but the probability can also be cal- 
culated by the formula: 


365 

Can you see why this formula produces the same result as the description 
in words and associated progression above? 

Consider all the presidents of the United States. (How many have there 
been to date?) Two of them, James Polk and Warren Harding, were bom on 
November 2. Is this to be expected given the size of the group? 

Since birthdays can be predicted, at least statistically, it ought to be pos- 
sible to predict deaths as well. It is interesting too that John Adams, James 
Monroe, and Thomas Jefferson all died on July 4th. Millard Fillmore and 
William Taft both died on March 8. What is the probability of that set of 
events? 
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10 PRINT “People Sate birthday probability" 
20 Q=364/365 
30 FOR N=2 TO 40 
40 P=100*( 1 — Q ) 

50 PRINT N, USING"##. ##" j P; : PRINT " X" 

60 Q=Q*(365-N) / 3 65 
70 NEXT N 


People Saie birthday probability 


2 

0.27 

X 

3 

0.82 

X 

4 

1.64 

X 

5 

2.71 

% 

6 

4. 05 

% 

7 

5.62 

X 

6 

7.43 

X 

9 

9.46 

% 

10 

11. 69 

X 

11 

14. 11 

X 

12 

16.70 

X 

13 

19.44 

X 

14 

22.31 

X 

15 

25.29 

X 

16 

28.36 

X 

17 

31.50 

X 

18 

34.69 

X 

19 

37.91 

X 

20 

41. 14 

X 

21 

44.37 

X 


22 

47.57 

X 

23 

58.73 

X 

24 

53.83 

X 

25 

56.87 

X 

26 

59.82 

X 

27 

62.69 

X 

28 

65.45 

X 

29 

68. 10 

X 

38 

70.63 

X 

31 

73.05 

X 

32 

75.33 

X 

33 

77.50 

X 

34 

79.53 

X 

35 

81.44 

X 

36 

83.22 

X 

37 

84.87 

X 

38 

86.41 

X 

39 

87.82 

X 

40 

89. 12 

X 


Now it is your turn to write a program. Here is a game to be played among 
your friends. You make a bet that they have similar preferences in colors. 
Each person in the group puts up one penny and you, because you are so sure 
of yourself (and probability), put up an amount in cents equal to the number 
of people in the group. If you win, you get to keep all of the money; if your 
friends win, each of them gets double his original bet, or two cents. 

Then each person selects a color (from a larger list than there are people, of 
course). If any two have picked the same color, you win; if all have picked 
different colors, they win. 

Since you would like to win, you need to know how many colors should be 
on the list for different size groups to give you a better than 50-50 chance of 
winning. You might like the probability to be up around 70% or so. You can 
produce the necessary table with a five-line program. Go to it! 
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Coins in a Pocket 


The next two programs were originally written by Glenda Lappan and 
M J. Winter and appeared in Creative Computing magazine. They are marvel- 
ous simulations for illustrating various aspects of probability. 

Coins in a Pocket is a simulation of the following situation. A newspaper 
costs 5 cents. A customer has 5 pennies and a dime in his pocket and offers to 
pay for his paper by letting you, the vendor, select at random two of the six 
coins. If you and this customer repeat this procedure for the 20 working days 
of a month, how much more or less than $1.00 (20 days x 5 cents) are you 
likely to have collected? 

The program below solves this little problem. The first random coin is se- 
lected from a group of six (Line 50 in which X = INT(RND*6)). The value 
of X can be 0, 1, 2, 3, 4, or 5. If it is 0, we assume the dime was selected and a 
second pick is not made, because it will surely be a penny. Thus 1 1 cents is 
added to the running total in Line 110. 

If the first coin is a penny (X = 1, 2, 3, 4, or 5) then we make a second pick 
from the five remaining coins. Again, if the dime is chosen (Y = 0), 1 1 cents 
is added to the total; otherwise 2 cents is added. 

As you can see from the sample runs, after a great number of trials, the 
average amount of money collected each day seems to be close to 5 cents. If 
you would like to convince yourself that the answer is, in fact, 5 cents, con- 
sider the following. Assign a letter to each coin, A - F. Let A be the dime and 
B through F be the pennies. Since all combinations are equally probable, here 
are all the possible combinations: 

AB 

AC BC 
AD BD CD 
AE BE CE DE 
AF BF CF DF EF 

There are five combinations with a dime (5x11 cents) and ten combina- 
tions of only pennies (10 x 2 cents). Add it up and you have 55 cents plus 20 
cents divided by a total of 1 5 combinations which equals an average value of 
75/15 or 5 cents. 
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10 PRINT "Simulates taking 
PRINT "from a dime and 
20 INPUT ''Number of trials" 
30 U = 0 


two coins at 
five pennies" 
j N 


random" 


40 FOR K=1 TO N 
50 X = INT ( R NO ( 1 )*6) 

60 IF X=0 THEN 110 
70 Y=INT( RND( 1 ) *6 ) 

80 IF Y=0 THEN 110 
90 U=U+2 
100 GOTO 120 
110 U=U+11 
120 NEXT K 

130 PRINT "Average value i s "; U/N; "cents " 


Simulates taking two coins at random 
from a dime and five pennies 
Number of trials? 10 
Average value Is 3.8 cents 


Simulates taking two coins at random 
from a dime and five pennies 
Number of trials? 100 
Average value Is 4.97 cents 


Simulates taking two coins' at random 
from a dime and five pennies 
Number of trials? 1000 
Average value Is 4.655 cents 
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Baseball Cards 


In statistics and probability, it is frequently necessary to predict the num- 
ber of trials on average until one is successful. For example, if you are trying 
to roll a four with one die, on average how many tries will it take until you do 
so? 

Since a die has six sides, the probability of rolling any number between 1 
and 6 is p = %. Thus, on any given roll, you have a % chance of rolling a 
four and % chance of rolling something else. This failure is designated q. So 
to be successful in rolling a four, we have a % chance on the first roll. On two 
rolls, the probability is two trials x the probability of one failure, then success 
(2 x % x l /e). On the third roll, the probability of success is three trials x the 
probability of two failures followed by one success (3 x % x % x Vs). 

Continuing this reasoning leads to the formula for the expected number of 
trials to success: 

E=lp + 2qxp + 3qxp + 4qx.p + ... 

This series can be reduced and solved for E which leads to: 


E = 


i-q 



So now the answer to the original problem can be calculated. The expected 
number of trials to roll a four is 1/p = 6 tries until success. 

But there is another way to approach this type of problem with the assis- 
tance of the computer. Consider the problem. Assume there are ten different 
prizes in Crunchies cereal boxes. How many boxes of cereal would you have 
to buy to obtain the complete set? The formula above can be expanded to 
solve this problem: 


N , N 
N-l 



N 


1 


+ -y = N(l+-^- + 


+- 


N 


For a value of ten, you can solve this formula by hand. However, let’s say 
you would like to solve this problem for baseball cards as found in packs of 
bubble gum. If there are 50 cards in a complete set, how many packs of gum 
must you buy, on average, to get a complete set. Write a computer program 
using the general formula above to solve for a set of any size. The answer for 
50 cards is 224.96 packs of gum; how many would you have to buy for a set of 
100 cards? 

There is another way of approaching this problem. This method is similar 
to that used to randomly select coins out of the pocket. In this case, the ran- 
dom group is set equal to the total number of cards; this value is accepted as 
input in Statement 20. Each purchase is set equal to a random number be- 
tween 1 and N (Statement 120). This card is then put into its proper place in 
the collection (Statement 130). However, if there is already a card there from 
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a previous purchase, we simply increment the purchase counter (Statement 
160) but do not get any closer to obtaining a complete set. After each pur- 
chase, we test to see if the set is complete (Statement 170), otherwise we go on 
buying more packs of bubble gum. 

When the set is complete, the number of packs of gum are tallied up and 
printed out. After a set of trials, the average is computed. This averaging 
value should be reasonably close to the value obtained by the formula al- 
though it will take a great number of trials before the two numbers are within 
1% of each other. 


10 PRINT "Simulates buying bubble gum": 
PRINT "with baseball cards Inside." 
20 INPUT "Cards In ser1es"jN 
30 INPUT "Nuiber of simulated trlals'jK 
40 DIM C ( lee) 

50 PRINT "Trial Packs" 

60 S=0 

70 FOR 1=1 TO K 
80 D = 0 

90 FOR J=1 TO N 
100 C ( 0 ) =0 
110 NEXT 0 

120 X*INT(RND(1)*N)+1 
130 C(X)=C(X)+1 
140 IF C ( X ) = 1 THEN 160 
150 GOTO 120 
166 D=D+1 

170 IF D=N THEN 190 
180 GOTO 120 
190 T=0 

200 FOR J=1 TO N 
210 T=T+C(0) 

220 NEXT J 
230 PRINT I , T 
240 S=S+T 
250 NEXT I 

268 PRINT "Average", S/K 


Simulates buying bubble gum 
with baseball cards Inside. 
Cards In series? 10 
Number of simulated trials 5 
Trial Packs 

1 35 

2 73 

3 29 

4 21 

5 24 

Average 36.4 
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Simulates buying bubble gum 
with baseball cards Inside. 
Cards In series? 50 
Number of simulated trials 5 
Trial Packs 


1 155 

2 167 

3 247 

4 212 

5 309 

Average 218 


If you run this second program for a set of 100 cards, it will sometimes run 
for a very long time before arriving at the answer. Try it with a set of 500 
cards as are used by some real bubble gum manufacturers. You can go have 
your dinner while the program computes just the first value. Be sure to 
change the dimension statement in Line 40 to C(500). 
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System Reliability 


As more and more people in the world come to depend upon mechanical 
and electronic devices in a myriad of different ways, it is important that these 
technological devices continue to function. 

Some years ago, the military came up with a measure that could be applied 
to all kinds of systems, big and small, to measure reliability. It is called 
“mean time between failures.” What this means is the length of time, on av- 
erage, between breakdowns. For a tank, this may be 100 hours, while for a 
spacecraft, the MTBF must be considerably longer than the planned mission. 

As we saw in an earlier section, it is frequently desirable to break down a 
large problem into smaller subproblems. To calculate the MTBF for a space- 
craft would be quite impossible. Instead, it is necessary to start with smaller 
systems and build up to the whole. 

Consider one of the electrical subsystems of a spacecraft. It uses five 
components as shown below, two parallel systems A and B, arranged in se- 
ries. The parallel subsystems are said to be “redundant.” This is one method 
of increasing reliability since the system will continue to work if at least one 
of the components works. 



A B 

If the manufacturer of the components stated that each one has a 60% 
probability of lasting 1000 hours, what is the chance that the entire system 
will last 1000 hours? Take a guess before reading on to the solution below. Is 
your guess greater than 60% or less? Why? 

The program below is a simulation of this system. Remember in subsystem 
A, it will continue to work if any of the three components works and in B if 
either of the two components works. However, the system will fail if all three 
of the A components or both of the B components fail. 

In Statement 30, the program is set up to make 500 trials of the system. 
Statement 50 selects a random integer between 1 and 10 for each component 
for each trial. If the value of this integer is 6 or under, the component is okay 
(60% probability of working at the end of 100 hours). If it is 7, 8, 9, or 10, 
this means the component has failed, and the program prints “Bad!” 
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10 PRINT "Simulates use of a circuit": 
PRINT “for 100 hours.": 

PRINT 

2 0 g- 0 

30 FOR 1=1 TO 500 

40 FOR J=1 TO 5 

50 C( J)=INT(RND(1)*10)+1 

60 IF C(0)>6 THEN 90 

70 C(J)=1 

60 60T0 100 

90 C(J)=0 

100 NEXT J 

110 IF C ( 1 ) =1 OR C ( 2 ) =1 OR C( 3 ) =1 THEN 140 
120 PRINT “Bad! "; 

130 GOTO 180 

140 IF C (4 ) = 1 OR C(5) = l THEN 170 
150 PRINT "Bad! 

166 GOTO 180 
170 S=S+1 
180 NEXT I 

190 R=INT(10000*(S/500)+.5)/100 
200 PRINT: 

PRINT “Reliability approx. "; R j " X* 

Simulates use of a circuit 
for 100 hours. 

Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 
Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! Bad! 

Reliability approx. 77.4 t 


Statement 1 10 tests whether subsystem A works and Statement 140 checks 
out subsystem B. If both work, a success is counted in Statement 170. 

There are, of course, shorter ways to do this. For example, all the individ- 
ual tests (Statements 110 to 140) can be replaced by one overall test (new 
Statement 110). 

110 IF (C(1)+C(2)+C(3))*(C(4)+C(5))=0 THEN 140 
130 S=S+1 
140 NEXT I 

190 R=INT(10000*{S/500)+.5)/100 


Simulates use of a circuit 
for 100 hours. 

Reliability approx. 79.4 t 

Now it is your turn to use what you have learned in this chapter and write 
a program to determine the mean time between failures of the system above. 
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Geometry and Calculus 


In this chapter several of the problem solving approaches from earlier 
chapters are used to solve geometric problems. You will find that many prob- 
lems can be solved in a variety of different ways — by applying a formula, by 
trial and error, by successive approximations, and, in some cases, by common 
sense. Perhaps the most important thing to learn from these problems and 
programs is how to analyze a problem to reach the solution quickly and 
painlessly. 


Crossed and Slipping Ladders 

Here is a simple problem of a slipping ladder. A ladder 25 feet long is 
placed so its foot is 7 feet from the base of a building. The base of the ladder 
slipped on some loose gravel so that the top is 4 feet lower than where it was 
to start. How far did the foot of the ladder slip? 



7 Y 


The diagrams show the two positions of the ladder. By the Pythagorean 
theorem we know that a 2 + b 2 = c 2 , hence, the equations needed to solve the 
problem are: 

x = V25 2 - 7 2 
y = v / 25 2 - (x-4) 2 

and the amount of slippage of the base is z = y - 7. These three equations are 
put into a computer program which quickly calculates an answer of 8 feet. 

10 PRINT "Slipping laddder" Slipping laddder 

20 C=25 Base slipped 8 ft. 

30 B = 7 

40 X=SQR(C*C-B*B) 

50 Y=SQR(C*C-(X-4)*(X-4) ) 

60 Z=Y-B 

70 PRINT "Base si 1 pped"Z"ft. " 

While the aritmetic in the above problem is not particularly messy, it is still 
no great joy to solve by hand. However, the computer is just as happy to do 
the problem with really messy dimensions, say a ladder length of 27.83 feet 
and a distance from the wall of 7.62 feet. 

20 C=27. 83 Slipping laddder 

30 B=7. 62 Base slipped 8.38614 ft. 
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Let’s consider an old problem found in many classic collections. Two lad- 
ders, one 25 feet long and the other 35 feet long lean against buildings on op- 
posite sides of an alley as shown below. The point at which the ladders cross 
is 12 feet above the ground. How wide is the alley? 



By using similar triangles twice we find that 

12/a + 12/b = 1 

Then, by applying the Pythagorean theorem and reducing, we obtain: 

a 2 - b 2 = 600 

Using one of the methods described in the Problem Solving chapter, you 
can solve these two simultaneous equations. Or they can be combined into 
one equation in which z, the width of the alley is 

z = a/ 35 2 - a 2 

By eliminating b, the following fouth degree equation is obtained 
a 4 - 24a 3 - 600a 2 4- 14400a - 86400 = 0 

Again, this can be solved using one of the methods in the chapter on Prob- 
lem Solving. This is the traditional approach, but there is another. The solu- 
tion is the intersection of the curves described by the two original equations. 
By using successive approximations, say in steps of 1 for a, you could solve 
for b t and b 2 in the following restated original equations 



12a 
a - 12 


b 2 = a/ 600 - a 2 

When b x falls below b 2 , reduce the step to 0. 1 to obtain a closer approxima- 
tion. By continuing this method of successive approximations, it is possible to 
obtain a very accurate solution. 


97 


Is there another approach? Y es, there is and it also avoids the quartic equa- 
tion. It uses the original equations in a trial and error procedure as described 
in the chapter on Sets and Repetitive Trials. See if you can write a program 
using this approach. 

Here is another classic problem for you to solve in any way that you like. 
What is the longest ladder that can be carried in a horizontal position around 
the corner made where a 1 2-foot wide alley meets one that is 8 feet wide. The 
diagram shows the problem. 
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Distance Between Coordinate Points 


This program solves for the distance between any two points in three- 
dimensional space defined by their x, y, and z coordinates. 

The formula for the distance between two points in three-dimensional 


space is: 

d = (x 2 - x^ 2 + (y 2 - yj 2 + (z 2 - z^ 2 

It would be desirable for the program to be able to solve any problem of 
this kind. One approach would be to use input statements to accept the point 
coordinates. Another is to use a data statement to define the points. Then 
only this one statement has to be changed for a new problem or set of 
problems. 

The program is written to calculate the distance between three sets of 

points. The points used were: 

0, 0, 0 and 3, 4, 5 

3.5, -4.7, 6.2 and -0.9, 3.0, 4.4 

67, 36, 82 and 54, 25, 90 

This calculation is used extensively in aerospace navigation. Can you deter- 
mine the angle or “compass heading’’ of the resulting flight path? Better start 
with just two dimensions. 


10 PRINT "Calculates distance between' 
20 PRINT "points In 3-0 space." 

30 PRINT: PRINT "Coordinates 
40 READ A, B, C. X, Y, 2 
50 L*SQR((A-X)"2+(B-Y)*2+(C-Z)"2) 

60 PRINT A; B; C 

70 PRINT X; Y ■, Z; TAB( 25) ; L 

80 PRINT "Hit any key to continue" 

90 A$=INKEY$: IF A$="" THEN 90 

100 GOTO 30 

110 DATA 0,0,0, 3,4, 5 

120 DATA 3. 5, -4. 7, 6. 2, -.9, 3, 4. 4 

130 DATA 67,36,82,54,25,90 


D1 stance" 


Calculates distance between 
points In 3-D space. 

Coordinates Distance 


0 0 0 
3 4 5 

Hit any key to continue 

Coordl nates 
3.5 -4.7 6.2 

-.9 3 4.4 

Hit any key to continue 

Coordl nates 
67 36 82 

54 25 90 

Hit any key to continue 


7.07107 

01 stance 
9.04931 

D1 stance 
18.8149 
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Area — by Calculation 

It is a simple matter to calculate the area of regular geometric figures by 
using the usual formulae. However, the problem becomes more difficult when 
it is necessary to calculate the area of two combined regular shapes or the 
area of irregular shapes. This program shows the method of analysis for solv- 
ing a problem of the first type, while the next program demonstrates four 
methods of dealing with irregular areas. 

The problem is to solve for the shaded area of the figure below for any 
value of the radius, R. 



The first step is to recall the formulae for calculating the area of a circle 
and square: 

A (circle) = pi * R 2 
A (square) = Side 2 or (2 * R) 2 

The difference in area between a square and a circle inscribed within its 
borders is: 

A (difference) = A (square) - A (circle) 
and the area of one corner is the difference divided by 4. The program below 
will calculate the area for any radius. 

Now it is your turn. Write a program to calculate the difference in area 
between a circle and square in which the square is inscribed within the circle. 
Now, change your program to calculate the difference in area for a triangle, a 
hexagon and an octagon inscribed within a circle. 
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10 PRINT "Area between cl rcl e" : PRINT "Inscribed In a square. 

20 INPUT "Radi us " ; R 

30 C=3. 14159*R*2 

40 S= ( 2*R ) *2 

50 D=S-C 

60 Dl=D/4 

70 PRINT "Trapped area*;D 
80 PRINT "One corner ";D1 


Area between circle 
Inscribed In a square. 
Radius? 10 
Trapped area 85.841 
One corner 21.4683 


Area between circle 
Inscribed In a square. 
Radius? 1 

Trapped area .85841 
One corner .214603 


Extend your program to calculate the difference in area for any regular fig- 
ure inscribed within a circle of radius 1.0. Set up a table or results as follows: 
Number of sides Difference in Area 

3 


What can you conclude from these results? Do these areas follow in some 
sort of progression? 


Area — by Integration 

In many cases it is necessary to determine the area of an irregular figure or 
the area under a curve where an exact formula is not available. The approach 
most commonly used is to divide up the enclosed area into small regularly 
shaped pieces and sum up the areas of all of these pieces. 

The easiest shape to use in these calculations is a rectangle. A group of 
rectangles can either be inscribed within the irregular figure or curve, or 
circumscribed around it. The first two diagrams show these two methods be- 
ing used to find the area under a curve. A third method is to use trapezoids 
which, depending upon the direction of curvature, will either be inscribed or 
circumscribed automatically. 
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A fourth method, known as Simpson’s Rule, essentially fits a senes of 
parabolas between the points of the curve and calculates the average area. It 
requires that the area be divided into an even number of parallel slits. Let us 
call the total number of divisions 2m which are h distance apart. The point on 
the curve where the first line intersects is y 0 , the second y,, and so on until 
y 2m . The area is then given by the formula: 

A= %h[(y 0 + y 2 m) + 4(yi + y3 + ■ • • 

+ y 2 m-i) + 2( yz + y* + y^- 2 > ] 

The area as calculated by Simpson’s Rule converges extremely quickly 
compared to the other methods. Nevertheless, as the number of intervals in- 
creases, all the methods approach the same limit. As might be expected, the 
trapezoidal approach converges more quickly than either of the approaches 
using rectangles. Compute the average of the two methods using rectangles. 

What do you get? Does this suggest another method? 

The methods used in this program involve the calculus. And you thought 

the calculus was difficult! Now you know otherwise. 


10 SCREEN 0,0 j ^ „ 

20 PRINT "Computes area under curve by^ 

30 PRINT "four methods of 1 ntegratlon. # 

40 PRINT "Define function In Line 480. # 

50 PRINT "Enter start and end values of x 
60 INPUT " (smaller one f 1 rst ) " ; A.B 
70 PRINT "Number Inscrlb Clrcum Trape- ^ra 
80 PRINT "Intvls Rctngls Rctngls zolds bolas 

90 M=-2 
100 S = 6 
110 M=M+3 

120 FOR N=M TO M+2 
130 C=S 
140 Q=0 

150 P=0 . 

160 D=2*N : ' Interval s 

170 H=(B-A)/D 
180 FOR 1=0 TO D-l 
190 X=A+I*H 

200 GOSUB 480 , „ , . , 

210 P=P+H*FN Inscribed Rectangle 

220 X=X+H 

240 Q=Q+H*FN 0 s' Circumscribed Rectangle 

258 NEXT I t . . 

260 T=(P+Q)/2 : ' Trapezoid 

270 X=A: GOSUB 480 
280 F1=FN 
290 X=B: GOSUB 480 
300 F2=FN 
310 U=F1+F2 

320 FOR J=2 TO ( 0-2 ) STEP 2 
330 X=A+0*H 
340 GOSUB 480 


Inscribed Rectartgle 

Circumscribed Rectangle 
Trapezol d 
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350 U=U+2*FN 
360 NEXT 0 
370 V-0 

380 FOR K=1 TO D-l STEP 2 

390 X=A+K*H 

400 GOSUB 480 

410 V=V+4*FN 

420 NEXT K 

430 S=(U+V)*(H/3) Parabola 

450 NEXT T N DjTAB(4,if>;TAB(13>;QjTAB(22):T;TAB ( 31 )i S 
460 IF D<64 THEN 110 

480 fJ, = {^ ((C " S)/( < C+S ) /2 )) >1E - 04 THEN 110 ELSE END 
490 RETURN 


Conputes area under curve by 
four lethods of Integration. 
Define function In Line 480. 
Enter start and end values of x 
(sialler one first)? 0 l 
Number Inscrlb Clrcui Trape- 
Intvls Rctngls Rctngls zolds 


.0625 
. 140625 
. 191406 
.219727 
. 234619 
.242249 


.5625 
. 390625 
.316406 
. 282227 
.265869 
.257874 


s zolds 
.3125 
.265625 
.253906 
.250977 
.250244 
. 250061 


Para- 
bolas 
. 583333 
.25 
.25 
.25 
. 25 
.25 


Cosputes area under curve by 
four iethods of Integration. 

Define function In Line 480. 

Enter start and end values of x 
(sialler one first)? 8 , 2 
Nuiber Inscrlb Clrcui Trape- Para- 
Intvls Rctngls Rctngls zolds bolas 


1 9 

2.25 6.25 

3.0625 5.0625 

3.51563 4.51563 

3.75391 4.25391 

3.87598 4.12598 


5 

4. 25 
4.0625 


4.51563 4.01563 
4.25391 4.00391 
4.12598 4.00098 


bol as 

9. 33333 

4 

4 

4 

4 

4 
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Science 


Using techniques and approaches presented in the previous chapters, this 
chapter contains five programs in the area of science. One uses a formula to 
solve simple gas problems, two are drill exercises on the gas laws of Boyle and 
Charles, and the last two are simulations. You will see that the simulations 
draw upon many previous techniques such as progressions and repetitive 
calculations. 


Gas Volume 


Here is a simple program to produce a not-so-simple table of values for gas 
volumes. 

The volume of a gas varies directly with the absolute temperature T (Kel- 
vin) and inversely with the pressure P. If a certain quantity of gas occupies 
500 cubic feet at a pressure of 53 pounds per square foot and an absolute tem- 
perature of 500 degrees, what volume will it occupy at 600 degrees absolute 
temperature and pressures from 100 to 1000 pounds per square foot in in- 
crements of 50 pounds? 

The original conditions are used to solve for the constant K in Statement 
50 (K = V*P/T). Then new volumes are computed for varying pressures 
with T equal to 600 degrees. The formula used is V = K*T/P. 

In the second part of the program, Lines 70 to 100 are replaced to produce 
a plot of the gas volume for the various pressures. 

How would you modify the program to deal with a more general case (i.e., 
other gasses and different temperatures)? Second, can you write a program 
that produces a table of values for a gas at different pressures and 
temperatures? 


10 PRINT “Gas volumes at different pressures." 

20 V=500 
30 P=53 
40 T=500 
50 K=V*P/T 

60 INPUT "Temperature (K)";T1 

70 PRINT "Pressure Volume Pressure Volume" 
80 FOR P=100 TO 1000 STEP 100 
90 V = K*T1 /P: V1 = K*T1 / ( P+50 ) 

100 PRINT P;TAB(10) } V ; TAB(20) i P+50 i TAB(30) i VI 
110 NEXT P 


Gas volumes at different pressures. 
Temperature (K)? 600 


Pressure 

Vol ume 

Pressure 

Vol ume 

100 

318 

150 

212 

200 

159 

250 

127.2 

300 

106 

350 

90.8571 

400 

79.5 

450 

70.6667 

500 

63.6 

550 

57.8182 

600 

53 

650 

48.9231 

700 

45.4286 

750 

42.4 

800 

39.75 

850 

37.4118 

900 

35.3333 

950 

33.4737 

1000 

31.8 

1050 

30.2857 
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10 PRINT "Gas volumes at different pressures." 

20 V=500 
30 P*53 
40 T=500 
50 K*V*P/T 

60 INPUT "Temperature (K)" ; T1 
70 CLS; SCREEN 0,0 

80 PRINT " Vol 0 200 400 600 800 1000"; 

90 PRINTiPRINT "220" : PRINT: PRINT "140": PRINT: PRINT "60" 

100 FOR P=100 TO 1000 STEP 10 

110 V*K*T1/P 

120 X=P/5+39 

130 Y*64-V/5 

140 IF Y<0 THEN Y=0 

150 PSET ( X, Y ) 

160 NEXT P 
170 GOTO 170 


Gas volumes at different pressures. 
Temperature (K)? 600 
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Charles' Law Drill 


In the previous program, gas volumes were calculated using Charles’ Law 
and Boyle’s Law. Here is a drill and practice program (remember Chapter 1?) 
that produces problems relating the volume and temperature of a gas. When 
pressure is constant, the volume/temperature relationships can be stated as 
follows: 

Vo = Vj 

To T, 

The program presents four problems, one to solve for each of the four vari- 
ables in the equation above. 

How can you make the program more efficient? More interesting? 


10 CLS: S = VAL(RIGHT$(TIME$, 2) ) 

20 PRINT TAB( 10 ) "Drl 1 1 on Charles' Law" 

30 PRINT" Vol umes are In ml and temp. In Kelvins." 

40 FOR D= 1 TO S 
50 V=INT(RND(1 )*50+50)*100 
60 V1=INT(RND{1)*50+50)*100 
70 T=INT(RND(1)*125+125) 

80 T1=INT(RND(1)*125+125) 

90 NEXT D 
100 C=C+1 

110 ON C GOTO 120,150,180,210,240 
120 V = 0-.PRINT "Solve for V when" 

130 Q1=V1*T/T1 
140 GOTO 250 

150 V1=0:PRINT “What Is VI given" 

160 Q1=V*T1/T 
170 GOTO 250 

180 T=0 : PRINT "Calculate the value of T" 

190 Q1=V*T1/V1 
200 GOTO 250 

210 T1=0:PRINT "Solve for T1 given" 

220 Q1=T*V1 / V 
230 GOTO 250 
240 STOP 

250 PRINT " V=" ; V; TAB( 10) " Vl = 11 ; VI; TAB(20 ) "T= " : T ; TAB( 30) "Tl = ": T1 

260 INPUT’Your Answer =";Q 

270 PRINT"The correct value=";Ql 

280 PRINT TAB(7)"Press any key to continue."; 

298 A$=INKEY$: IF A$>"" THEN CLS: GOTO 40 
300 GOTO 290 
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Drill on Charles' Law 
Volumes are In ml and temp. In Kelvins. 


Solve for V when 

V= 0 Vl= 9100 T= 151 Tl= 

Your Answer = ? 9300 

The correct value* 9476.55 

Press any key to continue. 

What Is VI given 

V= 7500 VI* 0 T= 202 Tl= 

Your Answer = ? 7600 

The correct value* 7648.52 

Press any key to continue. 

Calculate the value of T 

V* 7300 VI* 5400 T= 0 Tl= 

Your Answer = ? 280 

The correct value* 255.5 

Press any key to continue. 

Solve for T1 given 
V* 8400 VI* 6000 T= 140 Tl= 

Your Answer = ? 100 
The correct value* 100 

Press any key to continue. 

Break In 240 


145 


206 


189 


0 
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Boyle's Law Drill 


Boyle’s Law describes the behavior of gases under ideal conditions when 
pressure and volume are varied. When the temperature is constant, Boyle 
found that: 

P 0 * V 0 = Pi * V, 

Pressure is normally measured in centimeters of mercury while volume is 
measured in millilitres. As with the drill on Charles’ Law, this program 
presents four problems, one to solve for each of the four variables in the equa- 
tion above. 

Unlike many other drill and practice programs, these two do not compare 
your answer with the correct one. Instead, they leave that up to you to do. Is 
this desirable? Why or why not? If you feel it is undesirable, change the pro- 
grams so they do compare answers and calculate a score. 


10 CLS: S=VAL(RIGHT$(TIME$, 2) ) 

20 PRINT TAB( 9) "Drl 1 1 on Boyle's Law" 

30 PRINT "Voluraes are In ml and pressure in ram Hq" 

40 FOR 0=1 TO S 
50 V=INT(RND(1)*50+50)*100 
60 V1=INT(RND(1)*50+50)*100 
70 P=INT(RND(1)*150+150) 

80 P1=INT(RND(1)*150+150) 

90 NEXT 0 
100 C=C+1 

110 ON C GOTO 120,150,180,210,240 
120 V=0:PRINT "Solve for V when" 

130 Q1=P1*V1/P 
140 GOTO 250 

150 V1=0 : PRINT "What is VI given" 

160 Q1=P*V/P1 
170 GOTO 250 

180 P=0: PRINT "Calculate the value of P" 

190 Q1=P1*V1/V 
200 GOTO 250 

210 P1=0; PRINT "Solve for PI given" 

220 Q1=P*V/V1 
230 GOTO 250 
240 STOP 

250 PRINT "V= H i Vl;TA8(10)i " VI = " ; VI , TAB ( 20 ) ; "P= " ; P ; TAB ( 30 ) : "P1 = " ; P1 

260 INPUT"Your answer=";Q 

270 PRINT“Correct value="iQl 

280 PRINT TAB(7)"Press any key to continue."; 

290 A$ = INKEY$: IF A$>'"' THEN CLS:GOTO 40 
300 GOTO 290 
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Drill on Boyle's Law 
Volumes are in ml and pressure In mm Hg 


Solve for V when 
V= 8700 Vl= 8700 P= 
Your answer= ? 8700 
Correct value= 8793.21 
Press any key to 


280 PI* 283 


conti nue. 


What is VI given 
V= 0 Vl= 0 P* 

Your answer= ? 5600 
Correct value= 5696.24 
Press any key to 


163 PI* 186 
conti nue. 


Calculate the value 
V= 8500 Vl= 8500 
Your answer* ? 250 


of P 

P= 0 PI* 258 


Correct value* 284.805 

Press any key to continue. 


Solve for PI given 
V* 5300 VI* 5300 P= 155 PI* 0 
Your answer* ? 185 
Correct value* 187.17 

Press any key to continue. 


Break In 240 
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Photoelectric Emissions 

When light of a short wavelength falls on a metal surface, electrons are 
emitted from the metal. According to the description of this phenomenon by 
Einstein, there is a maximum wavelength for every metal above which no 
electrons are emitted. This is called the critical wavelength of the metal. 

This program simulates a laboratory experiment in which a metal is placed 
m a vacuum and bombarded with soft X-rays. The number of electrons emit- 
ted is collected and measured with a microammeter. The program simulates 
three trials at each of nine wave lenghts. 

After each set of experimental data, the program asks if you would like an- 
other run at a higher light intensity. The reason for doing this is that some- 
times at low light intensities, not enough electrons have been emitted for 
meaningful measurements. 

You can increase the precision of the experiment by increasing the number 
of wavelengths at which it is run. This value can be changed in Statement 70- 
note that the variable L is divided into 1000 in order to express the wave- 
length in Angstroms. One Angstrom (A) equals 10 centimeters or 10 
microns. 

Here are the coefficients for several metals: 


Silver .308 

Bismuth .338 

Cadmium .318 

Lead .340 

Platinum .385 


It is a rare physics laboratory in a high school or college today that has the 
experimental apparatus to run this experiment, yet with a small computer the 
equipment can be simulated. There are many other things with which you 
would not normally be able to experiment that can be simulated with a com- 
puter. Things such as a nuclear power plant, a malaria epidemic, an urban 
mass transit system, and a bicycle factory. 

Can you write a simulation for a real world system? You will find sections 
on simulations in the books Computers in Mathematics and Computers in Sci- 
ence and Social Studies . These, along with articles in Creative Computing 
magazine, might be of some help in writing a simulation of your own. 
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JJ P:?9: a " simulates metal bombardment with soft x-rays. 

20 PRINT: INPUT "Coefficient of metal = \V 
30 K=INT ( 1+2*RND( 1 ) ) 

40 CLS: PRINT "Output is In ml croamperes PRINT 
50 PRINT "Wave Trial Trial Trial" 

60 PRINT "Length 1 2 3"-PRINT 

70 FOR L= .42 TO .25 STEP-. 02 
80 M= INT( 1000/L ) 

90 PRINT M; 

100 FOR J=1 TO 3 

110 I=SQR(INT(25*RND(1))) 

120 IF L> V THEN 140 
130 GOTO 150 

140 I = SQR(K*K*100 + INT(35*RND(1) ) ) 

150 N=INT( 10*I+. 5 ) / 1 0 
160 PRINT TA8(J*8)iN: 

170 NEXT J 
180 PRINT 
190 NEXT L 

200 PRINT: INPUT"Increase light intensity (Y or N)" : AS 

210 IF A $ = " N " THEN 260 

220 IF A $ < > " Y " THEN 200 

230 INPUT"By what factor": F 

240 K=K*F 

250 PRINT: GOTO 50 

260 STOP 


This program simulates metal bombardment with soft x-rays. 


Coefficient of metal= ? .34 
Output Is In microamperes. 


Wave 

Trl al 

Trial 

Trial 

Length 

1 

2 

3 

2380 

20.2 

20.4 

20.7 

2500 

20.2 

20. 3 

20.7 

2631 

20. 1 

20.4 

20.3 

2777 

20 

20.6 

20.5 

2941 

4. 2 

1.7 

4.5 

3125 

0 

4.7 

4. 7 

3333 

3.9 

4.7 

4.2 

3571 

3.6 

4. 8 

3.3 

3846 

2.4 

3.2 

2.4 


Increase light Intensity (Y or N) ? Y 
By what factor ? 100 


Wave 

Trl al 

Trial 

Trial 

Length 

1 

2 

3 

2380 

2000 

2000 

2000 

2500 

2000 

2000 

2000 

2631 

2000 

2000 

2000 

2777 

2000 

2000 

2000 

2941 

3. 9 

2.4 

3.9 

3125 

2.2 

1.7 

3 

3333 

3.2 

1.7 

3.6 

3571 

2.8 

2.4 

4.5 

3846 

3 

3.5 

4. 5 


Increase light intensity (Y or N) ? N 
Break In 260 


Mutation of Moths 


This program is a simulation of the growth of a colony of pepper moths. 

The program allows a genetic mutation to be introduced in some year be- 
tween 1 and 30. The mutation can favor either dark or light colored moths. 

The program as it appears starts with a total colony of light moths; you 
could add an initial group of dark colored moths in Statement 45 as PI. 

The sample run shows a mutation which favors dark moths occurring in 
Year 3. The mutation affects about 2 % of the light moths each year and 
causes them to become dark. The program displays the number of moths of 
each color over a 30-year period. 

Obviously, because of the long time period involved it would be difficult to 
carry out this experiment in school, so the computer again is of real benefit. 

10 CLS:PRINT "This program simulates mutation of a":PRINT "colony of moth 
20 PRINT: INPUT "Rate of mutation (1 to 10)=";M 

30 PO=10000:Z = PO . L11 T41011TH H y 

40 PRINT “In what year does change in envi ronment INPUT occur ;X 

50 L- 1 : D=2 

60 PRINT "Change favor light or dark moths":INPUT "(L or D) M ; E$ 

70 PRINT 1 ' Year Dark Moths Light Moths" 

80 FOR T=1 TO 30 

90 IF T> = X THEN 120 

100 P1=0 

110 GOTO 170 

120 IF E $ < > 11 D “ THEN 100 

130 P1 = INT(P1 + . 01*M*PO+. 5) 

140 P0=INT (Z-P1+. 5) 

150 IF P1<Z THEN 170 
160 P1=Z:PO=0 

170 PRINT T; TAB( 10 ) PI t TAB ( 25 ) PO 
180 NEXT T 
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This program simulates mutation of a 
colony of moths. 

Rate of Mutation (1-10)= ? 2 

In what year does change In environment 

occur ? 3 

Change favor light or dark moths 


\ w ’ ' 
Year 

Dark Moths 

Light Moths 

1 

0 

10000 

2 

0 

10000 

3 

200 

9800 

4 

396 

9604 

5 

588 

9412 

6 

776 

9224 

7 

960 

9040 

6 

1141 

8859 

9 

1318 

8682 

10 

1492 

8508 

11 

1662 

8338 

12 

1829 

6171 

13 

1992 

8008 

14 

2152 

7848 

15 

2309 

7691 

16 

2463 

7537 

17 

2614 

7386 

18 

2762 

7238 

19 

2907 

7093 

20 

3049 

6951 

21 

3188 

6812 

22 

3324 

6676 

23 

3458 

6542 

24 

3589 

6411 

25 

3717 

6283 

26 

3843 

6157 

27 

3966 

6034 

28 

4087 

5913 

29 

4205 

5795 

30 

4321 

5679 


The so-called “killer bees” that were introduced into South America in the 
mid 70’s were supposed to help honey production because they were much 
more energetic than the lazy honey bees in Brazil. The idea was that they 
would mate with the existing honey bees and produce a more productive 
strain. However, they went on a rampage killing people and terrorizing the 
country. Then they started to migrate north. U.S. agricultural officials be- 
came alarmed that they would invade this country and bring death and 
destruction. 

Over the years, the killer bees have started to mate with the more docile 
South American honey bees, but only at the rate of 3% per year. Assuming 
they pose a danger to the U.S. until their numbers are reduced to 15% of 
their original quantity, how many years will it take for this to occur? The 
present migration patterns will bring them to the U.S. by 1989; will they still 
be dangerous (assume that the year of introduction was 1974)? 
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Projectile Motion 


The path followed by a projectile is called its trajectory. The trajectory is 
affected to a large extent by air resistance, which makes an exact analysis of 
the motion extremely complex. We shall, however, neglect the effects of air 
resistance and assume the motion takes place in empty space. 

In the general case of projectile motion, the body (bullet, rocket, mortar, 
etc.) is given an initial velocity at some angle 6 above (or below) the 
horizontal. 

If V 0 represents the initial velocity (muzzle velocity), the horizontal and 
vertical components are: 


V ox = V o cos0 > V oy = V o sin0 


Since we are neglecting air resistance, the horizontal velocity component 
remains constant throughout the motion. At any time, it is: 


V x = V ox = V o cos0=constant (0 

The vertical part of the motion is one of constant downward acceleration 
due to gravity. It is the same as for a body projected straight upward with an 
initial velocity V Q sin 6 . At a time “t” after the start, the vertical velocity is: 


V y ~ V oy ' = V o sin 6 - gt (2) 

where “g” is the acceleration due to gravity. 



(3) 


The horizontal distance is given by: 

x = V ox t = (V o cos0)t 

and the vertical distance by: 

y=V oy t-l/2 gt 2 

= (V o sin0)t- 1/2 gt 2 

The time for the projectile to return to its initial elevation is found from 
Equation (4) by setting y = 0. This gives 

2 V.sin 9 

t = 9 (5) 

g 

The horizontal distance when the projectile returns to its initial elevation is 
called the horizontal range. “R.” Introducing the time to reach the point in 
Equation (3), we find: 

2 V 2 sin 9 cos 6 

R = (6) 

g 

Since 2sin 9 cos 9 = sin 26, Equation 6 becomes: 

V 2 sin 29 

R = (7) 

g 

The horizontal range is thus proportional to the square of the initial veloc- 
ity for a given angle of elevation. Since the maximum value of sin 29 is 1, the 
maximum horizontal range, Rmax is V 0 2 / g. But if sin 29 = 1, then 29 90 

and 9 = 45°. Hence the maximum horizontal range, in the absence of air 
resistance, is attained with angle of elevation of 45°. 

From the standpoint of gunnery, what one usually wishes to know is what 
the angle of elevation should be for a given muzzle velocity v 0 in order to hit 
a target whose position is known. Assuming the target and gun are at the 
same elevation and the target is at a distance R, Equation (7) may be solved 
for 9. 

0 = 1/2 sin -1 

/ R \ (8) 

= 1/2Sin "(^) 

Provided R is less than the maximum range, this equation has two solu- 
tions for values of 6 between 0° and 90°. Either of the angles gives the same 
range. Of course, time of flight and maximum height reached are both greater 
for the high angle trajectory. 
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For example, say the maximum range of our gun is 10,000 yards and the 
target is at 5,900 yards: 


0 = 1/2 sin -1 


5,900 

10,000 


= 1/2 36° 


= 18°, or 90° - 18° = 72° 


Try the computer game Gunner which appears below. Use trial and error 
to try and destroy the target (see sample run). You get five chances per target. 

How many shots did it take to destroy all five targets? Did you ever fail to de- 
stroy a target in five trials? 

From the discussion above, you should realize that 45 degrees of elevation 
provides maximum range with values over or under 45° providing less range. 

The maximum range of the gun will vary between 20,000 and 60,000 yards 
and the burst radius of a shell is 100 yards. 

You can also determine the correct firing angle from sine tables, or a slide 
rule, a scientific calculator or a computer. You should be able to destroy ev- 
ery target with just one shot. What happens when the target is very close? 
Can you always use whole angles? 

Now write a computer program to accept the maximum range of your gun 
and the range to the target and then calculate the correct firing angle. You 
will have to solve two problems to write such a program: 

1. The Basic language does not have an ARCSIN function. However, the 
following formula may help. 

sin* 1 x = tan -1 ( — - — \ 

\/l^27 

2. You must convert from radians to degrees. 
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10 CLS: PRINT “You are commanding a gun crew. Hit" 

20 PRINT “within 100 yards of the target." 

40 GOSUB 530 

70 R=INT(40000*RND(1) +20000) 

80 PRINT "Max range" ; R; "yards" 

90 Z=0:S1=0 

110 T=INT(R*(.1+.8*RND{1))) 

120 S=0 : GOT 0 310 

140 PRINT "Min = 1 degree":G0T0 320 
160 PRINT "Max = 89 degrees GOTO 320 
180 PRINT "Over by ABS ( E) ;" yards ": GOTO 320 
200 PRINT "Short by ABS(E ); "yards "-. GOTO 320 
220 PRINT; PRINT "TARGET DESTROYED!" 

230 PRINT Si "rounds used";GOSUB 530 

250 S1=S1+S 

260 IF Z=4 THEN 450 

270 Z=Z+1 ; CLS 

290 PRINT "The forward observer has spotted More 
295 PRINT "enemy activity..." 

300 GOTO 110 

310 PRINT "Target at " ; T ( "yards " 

320 REM 

330 INPUT "El evatl on"; B 
340 If B> 89 THEN 160 
350 IF B<1 THEN 140 
360 S=S+1 

370 IF S<6 THEN 410 

380 CLS: PRINT TAB(20) "BOOM! PRINT TAB( 5) "BOOM! " 

385 PRINT TAB(27)"B00M! "-.PRINT TAB ( 10 ) "BOOM! GOSUB 530 
400 CLS:PRINT "The enemy got you first! ":GOTO 480 
410 B2 = . 035*B: I=R*SIN( B2 ) : X=T-I : E=INT( X ) 

420 IF ABS(E)<100 THEN 220 
430 IF E> 100 THEN 200 
440 GOTO 180 

450 CLS; PRINT "Total rounds = "jSl 

460 IF Sl> 18 THEN 480 

470 PRINT "Nice shootl ng! " : GOTO 490 

480 PRINT "Go back to Ft. Sill for more training." 

490 PRINT:INPUT "Try again (Y or N)" 5 Z$ 

500 IF Z$="y" OR Z$= " Y " THEN 70 
510 PRINT "Okay. Return to ca«p. " 

520 END 

530 FOR T= 1 TO 900:NEXT T 
540 RETURN 
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You are commanding a 
within 100 yards of t 
Max range 43808 yards 
Target at 8116 yards 
Elevation ? 10 
Over by 6906 yards 
Elevation ? 4 
Short by 2002 yards 
Elevation ? 5 
Short by 488 yards 
Elevation ? 6 
Over by 1017 yards 
Elevation ? 5.5 
Over by 266 yards 
Elevation ? 5.3 

BOOM! 

BOOM! 


The enemy got you first! 

Go back to Ft. Sill for lore training. 

Try again (Y or N ) ?Y 
Max range 50639 yards 
Target at 28461 yards 
Elevation ? 20 
Over by 4162 yards 
Elevation ? 19 
Over by 2787 yards 
Elevation ? 17 


gun crew. Hit 
ne target. 


TARGET DESTROYED! 

3 rounds used 

The forward observer has spotted more 
enemy actl vl ty . . . 

Target at 34829 yards 
Elevation ? 45 
Over by 15810 yards 
Elevation ? 35 
Over by 12813 yards 
Elevation ? 30 
Over by 9097 yards 
Elevation ? 25 
Over by 4039 yards 
Elevation ? 15 
Short by 9448 yards 
Elevation ? 17 

BOOM! 

BOOM! 


BOOM! 


BOOM! 


The enemy got you first! 

Go back to Ft. Sill for more tralnl 


ng. 


Try again (Y or N) ?N 
Okay. Return to camp. 
Break in 520 
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9 

Potpourri 


Here are five programs that didn’t seem to fit anywhere else. The first and 
the last are games, although you may come to think of the lunar landing 
simulation as a game also. One is a nifty simulation of smog, and the other 
calculates depreciation by three different methods. 
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Number Guessing Game 


Here is a computer program that plays the popular number guessing game. 
In it the computer picks a secret number between 1 and 100. You attempt to 
guess that number in as few tries as possible. 

There are many ways to go about guessing the secret number. Let some of 
your friends play this game and see what approaches they use to find the se- 
cret number. One approach is to start with a guess of 10. If this is too low, in- 
crease each guess by 10 until the computer says that a guess is too high. 
When this point is reached, start from the previous guess and increase each 
guess by 1. This is the method used to solve some of the problems in the chap- 
ter on convergence. 

Another approach is to try to bracket the number between upper and lower 
limits and reduce the limits by steps until the number is finally found. Two of 
the convergence programs used this approach. 

Is one of these the best way? Well, these methods are not bad, particularly 
compared to starting with 1 and simply counting to 100 until the solution is 
found. But there is a better way. It is known as binary search. 

This technique involves dividing the search domain, in this case 1 to 100, in 
half, and then in half again, and so on until the secret number is found. Play 
the game many times using different approaches. In the long run you should 
find that the binary search approach is the most efficient. 

In Line 230, the program contains the statement that “you should not need 
more than 7 guesses.” Why? If the secret number was between 1 and 128, 
what is the maximum number of guesses that would be necessary to find it? 
What if the number range were 1 to 130; then what would be the maximum 
number of guesses? 

Revise the program to choose a secret number between 1 and 10,000. Now 
the upper limit is 100 times the 1 to 100 game here which requires a maxi- 
mum of seven guesses to find the secret number; how many guesses will now 
be required? - 

Can you write a program in which the roles of the computer and player are 
reversed? In other words, the computer will try to guess your secret number 
between 1 and some upper limit. After each guess, you enter L for low, H for 
high, or C for correct. Can you write this program so the computer can tell if 
you are cheating, i.e., giving it inconsistent clues? 
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10 CLS: S=VAL(RIGHT$(TIME$, 2 ) ) 

20 PRINT "I've a secret number between 1 
25 PRINT "Try to quess It. I'll give you 
30 FOR 2=1 TO 1500: NEXT Z 
40 CLS 
50 C = 0 

55 FOR D=1 TO S 
60 N=INT(RND(1)*100)+1 
65 NEXT D 
70 C = C + 1 

80 PRINT TAB( 10 ) " Guess = " ; 

90 INPUT G 

100 IF N=G THEN 160 

110 IF G>N THEN 140 

120 PRINT "Too low.“:G0T0 70 

140 PRINT "Too high. “ : GOTO 70 

160 PRINT "Correct 1 n" ; Cj " tr 1 es . " 

170 IF 03 THEN 200 

180 PRINT "You were lucky. ":G0T0 240 

200 IF 07 THEN 230 

210 PRINT “Good Job.":G0T0 240 

230 PRINT "You shouldn't need more than 7 

240 FOR 2 = 1 TO 1500:.NEXT Z 

250 CLS 

260 INPUT “Play another game (Y or N)" ; A$ 
280 IF A $ = " Y * OR A$="y" THEN 20 
290 STOP 


I've a secret number between 1 and 100. 
Try to guess it. I'll give you clues. 
Guess= 50 

Too low. 

Guess= 75 

Too high. 

Guess= 62 

Too low. 

Guess= 67 

Too low. 

Guess= 73 

Too low. 

Guess= 74 
Correct In 6 tries. 

Good job. 

Play another game (Y or N) ? N 


Break In 290 


and 100. ■ 
cl ues . " 


tries. * 
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Depreciation — Three Methods 


This program shows how a piece of capital equipment depreciates accord- 
ing to three commonly used methods of depreciation: straight line, sum of the 
year digits, and double declining. 

The program asks for the original cost of the item, its expected life in years 
(the period of time over which it is to be depreciated), and its expected scrap 
(or sale) value at the end of that time. A table showing the annual depreci- 
ation for each of the three methods is then displayed. 


10 CLS:PRINT "Depreciation by three methods." 

20 LPRINT "Depreciation by three methods LPRINT " " 

30 INPUT "Original cost=";C 
40 LPRINT “Oriqinal cost=" ; C 
50 INPUT "Life in years="iL 
60 LPRINT "Life in years="jL 
70 INPUT "Scrap vaiue=";S 
80 LPRINT “Scrap Value= H ; S 
90 LPRINT ■ " 

100 V=C-S 
110 D1=V/L 

120 LPRINT “Straight line dep'n 1s ,, i Dl i "per year LPRINT" " 
130 Y-((L+1)/2)*L 
140 Z = L 

150 LPRINT " Sum of Double" 

160 LPRINT "Year Digits Declining" 

170 FOR X=1 TO L 
180 D2=V*( Z/Y ) 

190 D2=INT(D2*100+.5)/100 
200 Z=Z-1 
205 D3=2*C/L 

210 D3= INT ( D3*100+. 5) / 100 
230 C=C-D3 

240 LPRINT X ; TAB ( 5 ) D2 ? TAB ( 17 ) D3 
250 NEXT X 


Depreciation by three methods. 

Original cost= 8000 
Life in years= 8 
Scrap Value= 500 

Straight line dep'n is 937.5 per year. 



Sum of 

Doubl e 

Year 

Digits 

Decl 1 ni ng 

1 

1666. 67 

2000 

2 

1458. 33 

1500 

3 

1250 

1125 

4 

1041. 67 

843.75 

5 

833. 33 

632. 81 

6 

625 

474. 61 

7 

416.67 

355.96 

8 

208.33 

266.97 
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Smog Simulation 

This program is an adaptation of the smog model originally written by 
Herbert Peckham. The model assumes that vehicular traffic is the sole pro- 
ducer of smog, a somewhat poor assumption. It also assumes that most auto- 
mobile traffic occurs during the daylight hours and that traffic volume is very 
low (actually, zero) at night. The smog generated by the cars is dissipated by 
atmospheric conditions which vary depending upon sunlight, temperature, 
and weather. All of these conditions may be specified by the user. 

The model could be improved significantly by taking into account other 
sources of smog, by varying vehicular traffic according to the hour of the day, 
and by allowing daily variation of weather factors. Nevertheless, even in this 
rudimentary form it is interesting and instructive. 

A plot of the smog level is produced in Statement 410. Under some con- 
ditions, the smog level reaches a value that cannot be plotted because it is 
greater than the width of the screen (and printer). For runs with conditions 
like this, you might want to delete the plotting routine. A more elegant solu- 
tion would be to estimate the maximum value of the smog level from the in- 
put factors and calculate an appropriate plotting multiplication factor. 


10 CLS: 

PRINT "Smoq Simulation" 

20 PRINT 
30 T -7 

40 PRINT "Cars per road mile": 


PRIN1 "Los Angeles 200": 

PRINT "Detroit 100": 

PRINI "Tulsa 25" 

50 INPUT "Cars per m11e" ? C 
60 PRINT 

70 PRINT "Smog generation factor": 
PRINT "1950 Auto .005": 

PRINT "1975 Auto .001": 

PRINT "Bus .010" 


80 INPUT"Smog factorMl 
90 PRINT 

100 PRINT "Daytime smog breakdown": 
PRINT" 1 n percent per hour": 

PRINT “ Clear .02": 

PRINT " Cloudy .01" 

110 INPUT "Day breakdown"; R1 
120 PRINT 

130 PRINT "Nighttime breakdown": 
PRINT “ Hot .05": 

PRINT « Cool .10“ 

140 INPUT“Ni ght breakdown" ; R2 
150 PRINT 

160 PRINT "Dispersion (% per hour)": 
PRINT "High wind and rain 1.00' 
PRINT "No wind and dry .01* 
170 INPUT " D1 s per si on" ; R3 


180 PRINT 

190 PRINT “Table or plot": 

INPUT "Enter T or P" ; A$ 

200 PRINT: 

PRINT "Hour";TAB(22); "Smog Level 
210 R=R1 
220 K1=L1 

230 T 1 = INT( ( T-6 ) / 12 ) 

240 IF Tl/2=INT(Tl/2) THEN 270 
250 K1=0 
260 R=R2 

270 S=S+K1*C*10-R*S-R3*S 
280 IF S<=0 THEN S=0 
290 T=T+1 
300 T3=INT(T/12) 

310 T2=T-1 3*12+1 
320 X=INT(1000*S+.5)/1000 
330 IF T3 /2=INT ( T3/2 ) THEN 350 
340 PRINT USING"##" ; T 2 ; : 

PRINT" PM " ; T AB ( 10 ) ; : 

GOTO 360 

350 PRINT USING"##" ; T 2 ; : 

PRINT" AM" ; TAB( 10 ) ; : 

GOTO 360 

360 IF A$= “P " OR A$="p" THEN 390 
370 PRINT X 
380 GOTO 210 
390 X=4*X+5 

: 400 IF X> 38 THEN X=38 
410 PRINT TAB(X) 

420 GOTO 210 
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Smog Simulation 

Cars per road mile 
Los Anqeles 200 

Detroit 100 

Tulsa 25 

Cars per road mile ? 100 


Smog generation factor 


1950 Auto 

.005 







1975 Auto 

. 001 







Bus 

.010 







Smog factor ? 3E-03 


11 

AM 



* 




12 

AM 



« 

Daytime smog breakdown 


1 

PM 



* 

In percent 

per hour 


2 

PM 



* 

Cl ear 

.02 


3 

PM 



« 

Cl oudy 

.01 


4 

PM 



* 

Day breakdown ? .02 


5 

PM 



* 




6 

PM 



« 

Nighttime 

breakdown 


7 

PM 



♦ 

Hot 

. 05 


8 

PM 


« 


Cool 

. 10 


9 

PM 

* 



Night breakdown ? .05 


10 

PM 

• 






11 

PM 

* 



Di spersl on 

(% per hour) 


12 

PM 

* 



High wind 

and rain 1.00 


1 

AM 

* 



No wind and dry .01 


2 

AM 

* 



Di spersl on 

? .5 


3 

AM 

* 






4 

AM 

* 



Table or pi ot 


5 

AM 

* 



Enter T or 

P ? P 


6 

AM 

* 






7 

AM 

* 



Hour 

Smog 

Level 

8 

AM 


* 


9 AM 

* 


9 

AM 



* 

10 AM 

* 


10 

AM 



* 

11 AM 

# 


11 

AM 



• 

12 AM 


» 

12 

AM 



# 

1 PM 


» 

1 

PM 



* 

2 PM 


« 

2 

PM 



* 

3 PM 


« 

3 

PM 



* 

4 PM 


» 

4 

PM 



* 

5 PM 


« 

5 

PM 



♦ 

6 PM 


* 

6 

PM 



* 

7 PM 


* 

7 

PM 



* 

8 PM 

* 


8 

PM 


# 


9 PM 

* 


9 

PM 

* 



10 PM 

* 


10 

PM 

•* 



11 PM 

* 


11 

PM 

* 



12 PM 

# 


12 

PM 

* 



1 AM 

* 


1 

AM 

* 



2 AM 

« 


2 

AM 

* 



3 AM 

* 


3 

AM 

* 



4 AM 

* 


4 

AM 

* 



5 AM 

* 


5 

AM 

« 



6 AM 

* 


6 

AM 

* 



7 AM 

* 


7 

AM 

* 



8 AM 

* 


8 

AM 


* 


9 AM 

» 


9 

AM 



# 

10 AM 

•* 


10 

AM 



X 


127 


Smog City Simulation 


Cars per road mile 
Los Anqeles 200 

Detroit 100 

Tulsa 25 

Smog City 100 

Smog generation factor 
1950 Auto .005 

1975 Auto .001 

Bus .010 

Smog City IE-03 


Daytime smog breakdown 
in percent per hour 


Cl ear 

.02 


Cl oudy 

.01 


Smog City .01 

3 

Ni ghtti me 

breakdown 

4 

Hot 

. 05 

5 

Cool 

. 10 

6 

Smog City .1 

7 

8 

Dispersion (% per hour) 

9 

High wind 

and rai n 1. 00 

10 

No wind and dry .01 

11 

Smog City 

.05 

12 

1 



2 

Hour 

Smog Level 

3 

A 

9 AM 

1 

c 

10 AM 

1.94 

0 

11 AM 

2.824 

0 

7 

12 AM 

3.654 

o 

1 PM 

4.435 

o 

2 PM 

5. 169 

y 

i a 

3 PM 

5.859 

1 u 

1 i 

4 PM 

6. 507 

1 1 

5 PM 

7. 117 

12 

6 PM 

7.69 

1 

o 

7 PM 

8.228 

L 

8 PM 

6. 994 

J 

A 

9 PM 

5.945 

*» 

10 PM 

5.053 

b 

c 

11 PM 

4.295 

0 

*7 

12 PM 

3.651 

/ 

Q 

1 AM 

3. 103 

0 

2 AM 

2. 638 

y 

3 AM 

2.242 

10 

4 AM 

1.906 

1 1 

5 AM 

1. 62 

12 

6 AM 

1. 377 

1 

7 AM 

1.17 

2 

8 AM 

2. 1 

3 

9 AM 

2.974 

4 

10 AM 

3.796 

5 

11 AM 

4.568 

6 

12 AM 

5.294 

7 

1 PM 

5.976 

8 

2 PM 

6.618 

9 



10 


PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

PM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 

AM 


7.221 

7. 787 

8. 32 
8.821 

9. 292 
7.898 
6.713 
5.706 
4.85 
4. 123 
3.504 
2.979 

2. 532 
2.152 
1.829 
1. 555 
1.322 
2.242 

3. 108 
3.921 
4.686 
5.405 
6. 081 
6.716 
7.313 
7.874 
8.402 
8. 898 
9.364 
7.959 
6.765 
5.75 
4.888 

4. 155 
3. 532 
3.002 
2. 552 

2. 169 
1. 843 
1. 567 
1.332 
2.252 
3.117 

3. 93 
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Lunar Lander Simulation 

This program is one of the most popular computer simulations around. It 
is available in many versions; this one is adapted from the original program 
written in 1969. 

The program represents an exact simulation of an Apollo lunar landing 
module during the final descent. This portion of the descent would normally 
be controlled by the on-board computer backed up by another computer in 
the lunar orbiter, and still another computer on Earth. However, to exercise 
your knowledge of physics and to make an interesting game, we will assume 
that all three computers have had a simultaneous malfunction. Hence, it is up 
to you to land the spacecraft safely. 

To make a soft landing, you may change the bum rate of the retro rockets 
every ten seconds. You have a choice of not firing at all (bum rate of 0) or of 
firing at a fuel rate of between 8 and 200 pounds per second. Engine ignition 
occurs at 8 pounds, hence values between 1 and 7 pounds are not possible. 
You have 16,500 pounds of fuel. This is 500 pounds more than an actual 
LEM has, which will give you a little margin for error. When you get pro- 
ficient, change Statement 130 to N = 16000 to simulate the real thing more 
closely. The capsule weight is 33,000 pounds. 

Not that this is the way to come in, but if you did not fire the rockets at all, 
the estimated time for a free fall descent is 120 seconds to impact (and a huge 
splat). 

Good luck! 


10 CLS: 

PRINT TAB( 13) "LUNAR LANDER”: 

PRINT TAB ( 13 ) "by David Ah 1 * 

15 PRINT 

20 PRINT TAB(10) "Computer malfunction!": 

PRINT: 

PRINT" Take manual control of the capsule." 
30 G0SUB 810 
40 CLS: 

PRINT "Capsule weight Is 32,500 lbs.": 

PRINT "fuel weiqht Is 16,500 lbs." 

50 PRINT 

60 PRINT "Set retro burn rate every 10 seconds": 
PRINT “to zero or a value between 8 and 200": 
PRINT "lbs. per second." 

70 G0SUB 810 
80 PRINT 

90 PRINT "Time Miles M.P.H. Fuel Burn" 

95 PRINT STRINGS ( 38, ■»■) 

100 A=120 
110 V=1 
120 M=33000 
130 N=1 6500 
140 G=. 001 
150 2=1.8 
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160 P=INT( 3600*V+. 5) 

170 IF ABS( P ) > 1 THEN 190 
180 P=INT((3600*V+.5)*1000)/1000 
190 IF A<10 THEN 210 
200 D=INT(A+.5)i 
GOTO 240 

210 IF A< 1 THEN 230 
220 D=( INT( 10*A+. 5) ) / 10 
230 D=( INT( 1000*A+. 5) ) /1000 

240 PRINT INT(L); TAB(6) D; TAB(14) P, TAB ( 22 ) M-N; TAB(31): 

250 INPUT K 

260 IF K> 200 THEN 290 

270 IF K> 7 THEN 310 

280 IF K=0 THEN 310 

290 PRINT "Can't do - again please": 

GOTO 250 
310 T = 1 0 

320 IF M-N< . 001 THEN 430 
330 IF T < . 001 THEN 160 
340 S=T 

350 IF M> N+S*K THEN 370 
360 S= (M-N) /K 
370 GOSUB 770 
380 IF I<*0 THEN 630 
390 IF V<=0 THEN 410 
400 IF J<0 THEN 690 
410 GOSUB 570 
420 GOTO 320 

430 PRINT "You ran out of fuel at "; INT( L*1000 )/ 1000; "seconds. " 

440 S=(-V+SQR(V*V+2*A*G) )/G 

450 V=V+G*S 

460 L = L+S 

470 W=3600*V 

480 PRINT: 

PRINT “On noon at"; I NT ( L ) ; "seconds": 

PRINT "Impact velocity was • ; INT ( 100*W+. 5) / 100; "MPH. " 

490 IF W> 1 . 2 THEN 510 

500 PRINT "Perfect landing!"; 

END 

510 IF W>10 THEN 530 

520 PRINT "Oooooaph! That's rough on the shocks.": 

PRINT "You need more practice!": 

END 

530 IF W> 60 THEN 550 

540 PRINT "Severe craft damage!": 

PRINT "You're stranded until help arrives.": 

END 

550 PRINT "Sorry, there were no survivors.*: 

PRINT "You blew It. " 

560 PRINT "In fact, you blasted a new lunar": 

PRINT"crater"; INT(W*. 278); "feet deep! ": 

END 

570 L=L+S 
580 T=T-S 
590 M=M-S*K 
600 A=I 
610 V=J 
620 RETURN 

630 IF S< . 005 THEN 470 

640 D=V+SQR( V*V+2*A*(G-Z*K/M) ) 

650 S=2*A/D 

660 GOSUB 770 

670 GOSUB 570 

680 GOTO 630 

690 W=(l-M*G/(Z*K))/2 
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700 S=M*V/ (Z*K*(W+SQR(W*W+V/2) ) )♦. 05 

710 GOSUB 770 

720 IF I <=0 THEN 630 

730 GOSUB 570 

740 IF J>0 THEN 320 

750 IF V>0 THEN 690 

760 GOTO 320 

770 Q=S*K/M 

780 0=V+G*S-Z*Q*( 1+Q* ( . 5+Q* ( 1 /3+Q* ( . 25+Q/5) ) ) ) 

790 I=A-G*S*S/2-V*S+Z*S*Q*( . 5+Q*( l/6+Q*(l/12+Q/20) ) ) 

800 RETURN 

810 FOR H=1 TO 1250: LUNAR LANOER 

NEXT H: by Ahl 

RETURN 

Computer malfunction! 

Take manual control of the capsule. 


Capsule weight Is 32,500 lbs. 
fuel weight Is 16,500 lbs. 

Set retro burn rate every 10 seconds 
to zero or a value between 8 and 200 
lbs. per second. 


T 1 me 

Miles 

N.P.H. 

Fuel 

Burn 

0 

120 

3600 

16500 

? 0 

10 

110 

3636 

16500 

? 0 

20 

100 

3672 

16500 

? 25 

30 

90 

3659 

16250 

? 25 

40 

79 

3645 

16000 

? 25 

50 

69 

3631 

15750 

? 50 

60 

59 

3566 

15250 

? 100 

70 

50 

3394 

14250 

? 100 

80 

41 

3216 

13250 

? 100 

90 

32 

3031 

12250 

? 100 

100 

24 

2837 

11250 

? 200 

110 

16 

2388 

9250 

? 200 

120 

10 

1901 

7250 

? 200 

130 

5.909 

1367 

5250 

? 200 

140 

2.917 

777 

3250 

? 200 

150 

1. 652 

122 

1250 

? 200 

You ran out of 

fuel at 

156.25 

seconds 


On moon at 357 seconds 
Impact velocity was 394.53 MPH. 
Sorry, there were no survivors. 
You blew It. 
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Hammurabi 

Hammurabi is one of the all-time favorite computer games. On the one 
hand, it may be considered a game, but on the other it is an intriguing simula- 
tion of barter and management. 

Hammurabi is your servant as you try to manage the ancient city-state of 
Sumeria. The economy of the city-state revolves around just one thing — the 
annual crop of grain (probably soybeans). 

Each year, you must determine how many bushels of grain you wish to 
feed to your people (you’ll quickly discover how much a person needs to sur- 
vive), how much you wish to use as seed in planting crops for the coming 
year, how much you wish to use for the purchase of additional land from 
your neighboring city-state, and how much you wish to put in storage. 

Of course, if you have a bad harvest or if rats overrun your grain storage 
bins, you may have to sell land in order to get enough grain to keep your peo- 
ple from starving, or to plant the land for the coming year. Unfortunately, 
disasters always seem to strike, forcing you to sell land, when the price is at 
an all-time low; but that’s not any different from the real world. 

Most people start to play this game with noble ambitions. However, before 
long, they start longing for a plague to trim their growing population. Or they 
deliberately starve some people to keep things in balance (gosh, maybe these 
zero population growth people have something, after all!). 

Over the years, this game more than any other, has spawned a host of look- 
alikes, extensions, and modifications. Indeed, several manufacturers have 
taken my original with no changes whatsoever, put it in a fancy box, and 
charged a handsome price for it. Accept no imitations! Here is the original 
game (with the dialog shortened slightly) for you to run on your computer. 

If you want to experiment with changes, here are some suggestions. In the 
existing game, plagues randomly occur 15% of the time; lower this to 10% or 
5%. People now require a fixed amount of food; vary this amount slightly 
from year to year. Permit the construction of a rat-proof grain bin, but this 
must cost a fair amount. Introduce a mining industry as well as agriculture. 
How about fishing or tourism? Let your imagination run wild. Experiment! 
Have fun! 
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1 CLS: SCREEN 0,0 

2 PRINT TAB( 13) "Hammurabi Game" : PRINT 

10 PRINT "Try to govern ancient Sumerla for ten":PRINT 
20 FOR X=1 TO 1250: NEXT X 
30 01=0: P1=0 

40 2=0:P=95:S=2800:H=3000:E=H-S 
50 Y=3:A=H/Y: 1=5: Q=1 
60 D=0 

70 CLS:Z=Z+1 

80 PRINT "Hammurabi: Year";Z 
90 PRINT D| “starved, Ij "born; 

100 P = P+I 

110 IF Q>0 THEN 140 
120 P=INT(P/2) 

130 PRINT "Horrible plague! " 

140 PRINT "population =";P 

150 PRINT "You own" ; A; "acres and harvested"; Y 
160 PRINT "bushels/acre. Rats ate" ; E; "bushel s . * 

190 PRINT S; "bushels stored." 

200 IF Z=ll THEN 840 

210 C=INT(10*RND(1)):Y=C+17 

220 PRINT "Land Is trading at"; Y ; "bushels/acre. " 

240 INPUT'Acres to buy";Q 
250 IF Q<0 THEN 810 
260 IF Y*Q<S THEN 290 
270 GOSUB 730 
280 GOTO 240 
290 IF Q=0 THEN 330 
300 A=A+Q:S=S-Y*Q:C=0 
310 GOTO 390 

330 PRINT026O, “Acres to sell";: INPUT Q: IF Q<0 THEN 810 

340 IF Q<A THEN 370 

350 GOSUB 760 

360 GOTO 330 

370 A=A-Q: S=S+Y*Q: C=0 

390 INPUT"Bushel s for feed" ; Q 

400 IF Q<0 THEN 810 

410 IF Q<=S THEN 440 

420 GOSUB 730 

430 GOTO 390 

440 S=S-Q:C=1 

460 INPUT"Acres to seed";D:IF D=0 THEN 570 
470 IF D<0 THEN 810 
480 IF D<=A THEN 500 
490 GOTO 460 

500 IF INT ( D/2 ) <=S THEN 530 

510 GOSUB 730 

520 GOTO 460 

530 IF D<10*P THEN 560 

540 PRINT "Not enough people" 

550 GOTO 460 
560 S=S-INT (D/2) 

570 GOSUB 790 
580 Y=C: H=D*Y ; E=0 
590 GOSUB 790 

600 IF INT(C/2)<>C/2 THEN 620 
610 E=INT( S/C ) 

620 S=S-E+H 
630 GOSUB 790 

640 I=INT(C*(20*A+S) /P/100+1) 

650 C=INT(Q/20 ) 

660 Q=INT(10*(2*RND(1)-. 3) ) 


"years. " 
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670 IF P<C THIN 60 
680 D=P-C:IF D> . 45*P THEN 710 
690 PI = ( (2-l)*Pl+D*100/P)/Z 
700 P=C:D1=D1+D:00T0 70 
710 PRINT Di "people starved!" 

720 PRINT "You are 1 mpeached. " ; GOTO 1030 

730 PRINT “Think again. You have just"; S; "bushels" 

750 RETURN 

760 PRINT "Think aqaln. You own onl y " ? A; "acres 
780 RETURN 

790 C=INT (RND( 1 )*5)+l 
800 RETURN 

810 PRINT "Hammurabi : I can't do that. Get" 

820 PRINT “yourself another steward!" 

830 GOTO 1030 

840 PRINT "In ten years" ; Plj"% starved per year." 
850 PRINT "A total of “ : Dl ; * peopl e died!" 

870 L=A/P 

880 PRINT “You started with 10 acres/person and 
890 PRINT "ended wl th ", Li "acres /person. " 

910 PRINT 

920 IF Pl>33 THEN 720 
930 IF L<7 THEN 720 
940 IF Pl> 10 THEN 990 
950 IF L<9 THEN 990 
960 IF PI > 3 THEN 1010 
970 IF L<10 THEN 1010 

980 PRINT "Fantastic per for mance! GOTO 1030 
990 PRINT "Very heavy handed!" 

1000 PRINT "People are rebel 1 1 ng! GOTO 1030 
1010 PRINT "Not bad, but could be somewhat better." 
1030 PRINT "So long for now." 

1050 END 
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Hammurabi Game 

Try to govern ancient Sumerla for ten 
years . 


Hammurabi: Year 1 

0 starved, 5 born; population = 100 
You own 1000 acres and harvested 3 
bushels/acre. Rats ate 200 bushels. 
2800 bushels stored. 

Land Is trading at 22 bushels/acre. 
Acres to buy ? 0 Acres to sell ? 0 

Bushels for feed ? 2000 
Acres to seed ? 990 


Hammurabi: Year 2 

0 starved, 7 born; population = 107 
You own 1000 acres and harvested 1 
bushels/acre. Rats ate 76 bushels. 

1219 bushels stored. 

Land Is trading at 18 bushel s /acre. 

Acres to buy ? 0 Acres to sell ? 100 
Bushels for feed ? 2000 
Acres to seed ? 900 


Hammurabi: Year 3 

7 starved, 8 born; population = 108 
You own 900 acres and harvested 2 
bushels/acre. Rats ate 0 bushels. 

2369 bushels stored. 

Land Is trading at 25 bushel s /acre. 
Acres to buy ? 0 Acres to sell ? 0 

Bushels for feed ? 2000 
Acres to seed ? 900 

Think again. You have just 369 bushels 
Acres to seed ? 369 


Hammurabi : Year 4 

8 starved, 2 born; population * 102 
You own 900 acres and harvested 3 
bushels/acre. Rats ate 0 bushels. 

1292 bushels stored. 

Land Is trading at 24 bushel s /acre. 

Acres to buy ? 0 Acres to sell ? 50 

Bushels for feed ? 2000 
Acres to seed ? 850 


Hammurabi: Year 5 

2 starved, 8 born; population = 108 
You own 850 acres and harvested 1 
bushels/acre. Rats ate 0 bushels. 

917 bushels stored. 

Land Is trading at 25 bushels/acre. 

Acres to buy ? 0 Acres to sell ? 50 

Bushels for feed ? 1750 
Acres to seed ? 850 
Acres to seed ? 800 
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Hammurabi : Year 6 

21 starved, 2 born ; population = 89 
You own 800 acres and harvested 5 
bushels/acre. Rats ate 0 bushels. 

4017 bushels stored. 

Land Is trading at 20 bushel s /acre. 
Acres to buy ? 75 
Bushels for feed ? 1750 
Acres to seed ? 850 


Hammurabi : Year 7 

2 starved, 9 born; Horrible plague! 
population = 48 

You own 875 acres and harvested 1 
bushels/acre. Rats ate 171 bushels. 

1021 bushels stored. 

Land is trading at 19 bushel s /acre. 

Acres to buy ? 0 Acres to sell ? 75 

Bushels for feed ? 1000 

Acres to seed ? 480 

Not enough people 

Acres to seed ? 470 


Hammurabi : Year 8 

0 starved, 12 born; population = 60 
You own 800 acres and harvested 3 
bushel s /acre. Rats ate 0 bushels. 

2621 bushels stored. 

Land Is trading at 23 bushel s /acre. 
Acres, to buy ? 0 Acres to sell ? 0 

Bushels for feed ? 1100 
Acres to seed ? 590 


Hammurabi: Year 9 

5 starved, 3 born; Horrible plague! 
population = 29 

You own 800 acres and harvested 1 
bushels/acre. Rats ate 0 bushels. 

1816 bushels stored. 

Land Is trading at 18 bushel s /acre. 
Acres to buy ? 0 Acres to sel 1 ? 0 

Bushels for feed ? 600 
Acres to seed ? 280 


Hammurabi: Year 10 
0 starved, 19 born; population = 48 
You own 800 acres and harvested 2 
bushels/acre. Rats ate 0 bushels. 

1636 bushels stored. 

Land Is trading at 20 bushel s /acre. 
Acres to buy ? 0 Acres to sell ? 0 

Bushels for feed ? 975 
Acres to seed ? 470 
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Hammurabi : Year 11 

0 starved, 11 born; Horrible plague! 
population = 29 

You own 800 acres and harvested 2 
bushel s /acre. Rats ate 0 bushels. 

1366 bushels stored. 

In ten years 5.8727471767699 % starved per year. 
A total of 45 people died! 

You started with 10 acres/person and 
ended with 27.586206896552 acres/person. 

Not bad, but could be somewhat better. 

So long for now. 
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The magazine referred to in the text is: 

• Creative Computing. This is the leading magazine of software and applica- 
tions for all small computers. It carries articles, tutorials, how-to applications, 
and extensive in-depth evaluations. 
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• Computers in Science and Social Studies. Scores of simulation programs in 
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Write or call for the current price: 

Creative Computing 
39 E. Hanover Ave. 

Morris Plains, NJ 07950 
(800) 631-8112 
In NJ (201) 540-0445 
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THE TRS-80 MODEL TOO 



BY DAVID H. AHL 


i -- ■ ; 

This Ideabook contains dozens of ways to make the most out of your 
computer for solving practical, everyday problems. The 50 ready-to-run 
programs demonstrate scores of different techniques for solving problems in 
mathematics, science'and business. 


The ten chapters deal with solving problems by formulas and repetitive 
* t nal&^o nvergence, recursion, compounding, probability, geometry, science, 
simulations, and drill and practice. 


Some of the problems demonstrate the capabilities of the computer; 
others identify its shortcomings. It is important to be familiar with both the 
strengths and weaknesses of your tools so you can recognize the types of 
pbs for which they are suitable. ' ' • 


The author, David H. Ahl, has been involved with the use of com 
puters since 1957. fHe is the author of 16 books and is the founder and 
editorial director of Creative Computing, SYNC and Video & Arcade Games 


magazines. 



■ CREATIVE 

i COMPUTING 

■ PNFHS 


ISBN 0 916688-57-7 


